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The semileptonic decays of Ac and A(, are treated in the framework of a constituent quark model. 

Both nonrelativistic and semirelativistic Hamiltonians are used to obtain the baryon wave functions 
from a fit to the spectra, and the wave functions are expanded in both the harmonic oscillator and 
Sturmian bases. The latter basis leads to form factors in which the kinematic dependence on is in 
the form of multipoles, and the resulting form factors fall faster as a function of in the available 
kinematic ranges. As a result, decay rates obtained in the two models using the Sturmian basis 
are significantly smaller than those obtained using the harmonic oscillator basis. In the case of the 
Ac, decay rates calculated using the Sturmian basis are closer to the experimentally reported rates. 

However, we find a semileptonic branching fraction for the Ac to decay to excited A* states of 11% 
to 19%, in contradiction with what is assumed in available experimental analyses. Our prediction 
for the As semileptonic decays is that decays to the ground state Ac provide a little less than 70% 
of the total semileptonic decay rate. For the decays At, —> Ac, the analytic form factors we obtain 
satisfy the relations expected from heavy-quark effective theory at the non-recoil point, at leading 
and next-to-leading orders in the heavy-quark expansion. In addition, some features of the heavy- 
quark limit are shown to naturally persist as the mass of the heavy quark in the daughter baryon is 
decreased. 
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I. INTRODUCTION AND MOTIVATION 

Many of the parameters of the Standard Model (SM), including the Cabbibo-Kobayashi-Maskawa (CKM) P| 
matrix elements, are not yet determined with ‘satisfactory’ precision. Very precise knowledge of these matrix 
elements is important as they play a crucial role in the search for answers to some fundamental questions, such 
as the nature of CP violation and the unitarity of the CKM matrix. Semileptonic decays of hadrons have been, 
and will continue to be, the main source of information on the CKM matrix elements. The precision with which 
the CKM matrix elements are extracted from these semileptonic decays is strongly dependent on how well the 
form factors that describe the matrix elements of the hadronic currents are known. The vast literature on these 
form factors is a testament to the importance of these parameters. 

The semileptonic decays of heavy mesons have been studied extensively in the last two decades. Wirbel, 
Stech and Bauer Hi assumed monopole type form factors for the decays of heavy mesons. In Ref. i i, 
a non-relativistic quark model (NRQM) was used to treat the semileptonic decays of B and D mesons, and 
relatively simple forms for the form factors were presented. The first of those articles, along with the work of 
Shifman and Voloshin |^, ultimately lead to the development of the heavy quark effective theory (HQET). In 
addition, Ivanov and Santorelli 0 used a relativistic quark model to find the form factors. These are just a 
very few of the very large number of articles that treat semileptonic decays of mesons in some kind of model. 

Weak decays of hadrons involving one or more heavy quarks (mg >> Aqcd) have an additional symmetry 
in the effective Lagrangian which was first pointed out by Isgur and Wise [g. There, they used the additional 
heavy quark symmetry to obtain normalized, model-independent predictions for all the form factors for the 
decays of heavy hadrons to daughter hadrons that are also heavy. This led to many subsequent calculations by 
many authors. 
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For the hadronic matrix elements of the electroweak currents between two heavy mesons, the application of 
HQET provides a number of features that simplify the extraction of CKM matrix elements from such decays. 
First, the number of form factors is reduced, so that the six form factors that describe the decays of heavy pseu¬ 
doscalar mesons to heavy pseudoscalar and vector mesons are replaced by a single form factor, at leading order 
in HQET. This form factor has become known as the Isgur-Wise function. Second, the absolute normalization 
of this form factor at the so-called non-recoil point is known. Third, corrections to this normalization do not 
arise at order l/mg in the heavy quark expansion, but at order l/mg. This is known as Luke’s Theorem 
and is an analog of the Ademollo-Gatto theorem ^^1 ■ This means that some predictions made at leading order 
are more robust than might be expected. Finally, the corrections that do arise can be estimated systematical^ 
in the heavy quark expansion. As a result, HQET has become the tool of choice in the extraction of |14h| [ij- 

For the semileptonic decays of a heavy meson to a light meson, the predictions of HQET are not quite 
as powerful: there is no reduction in the number of form factors needed to describe the decay, nor are the 
normalizations of any of the form factors known. However, the heavy quark symmetry, along with SU(2) or 
SU(3) flavor symmetry for the light mesons, can be used to relate the form factors for H ^ tt and D ^ p decays, 
for instance, to those for H —> tt and B ^ p decays, respectively. Thus, even though it is not as predictive 
in the decays of heavy to light mesons, there is still a great deal of reliance on HQET for extracting Vub from 
meson decays. 

For the semileptonic decays of a heavy baryon to another heavy baryon, HQET makes predictions that are 
completely analogous to those made for heavy-to-heavy meson decays: (i) the six form factors that describe the 
decays to the ground-state heavy baryons are replaced by a single form factor, the Isgur-Wise function; (ii) the 
normalization of the Isgur-Wise function is known at the non-recoil point; (iii) corrections to this normalization 
first appear at order 1 /toq; (iv) corrections can be systematically estimated in a l/mg expansion. 

In the case of a heavy baryon decaying to a light baryon, HQET makes predictions that are not as powerful 
as in the heavy-to-heavy case, but which are significantly more powerful than for the heavy-to-light transitions 
of mesons. Among the baryons, the leading-order prediction is that the number of independent form factors 
decreases from six to two. In addition, as with mesons, the heavy quark symmetry can be used to relate the 
form factors for the decay A+ ^ n to those for A;, ^ p, for instance. This, in principle, could facilitate the 
extraction of Vub from semileptonic decays of the A},, and since the number of unknown form factors is reduced 
from six to two, the theoretical uncertainty in the extraction from these decays should be significantly smaller 
than extractions from meson decays. 

While HQET has been tremendously successful and useful in treating semileptonic decays of heavy hadrons, 
it is not without its limitations. It is a limit of QCD that applies only to hadrons containing heavy quarks. 
For the decays of such hadrons, it only predicts the relationships among form factors, not their kinematic 
dependence; ansatze, models of one kind or another, or lattice simulations, are still needed for this. In addition, 
the predictions of HQET are valid only as long as the energy of the daughter hadron is not comparable to the 
mass of the heavy quark. For heavy to heavy decays, this means that the predictions are valid for all of the 
available phase space, but for heavy to light decays, such as H ^ tt, a large portion of the available phase space 
is beyond the region of reliable applicability of HQET. These limitations mean that the predictions of HQET 
must be complemented/supplemented by information arising from other approaches to hadron structure. 

While some work has been done in modeling the form factors for the semileptonic decays of heavy baryons, 
to the best of our knowledge little has been done in treating the decays to excited baryons. Predictions for the 
number of independent form factors for decays to excited states have been made in the framework of HQET , 
and Leibovich and Stewart have examined the form factors for decays to the 1/2“ and 3/2“ states, using 
large Nc arguments. In the semileptonic decays of B mesons to those with charm, it is known that decays to 
the ground state pseudoscalar and vector mesons provide only about 75% of the total semileptonic decay rate, 
while for the B^, the corresponding fraction is about 85%. Any assumption that decay of a heavy baryon to 
the ground state will saturate the semileptonic decay rate is therefore subject to potentially large corrections. 

In some of the work done in this area, the predictions of HQET, along with various ansatze for the form 
factors, have been used to estimate some decay rates. Leibovich and Stewart ^ follow such a procedure to 
estimate the rates for decays of the A;, to the = 1/2“ and 3/2“ A^ states. Polarization effects in semileptonic 
Ab and Ac decays have been studied by Korner and Kramer Hi, using the predictions of HQET to estimate 
the dominant form factors for both 6 —> c and c —> s transitions. They have also calculated the asymmetry 
parameters that characterize the angular dependence of the decay distributions. 

A number of authors have constructed explicit quark models of the form factors for the decays of Ag baryons 
to ground state baryons. The decays of A^ and H;, have been treated by Singleton using a spectator quark 
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model 0 He also discusses the polarization of the W boson and the daughter baryon in these processes. 
Albertus et al. use a NRQM to evaluate the form factors for Ah ^ Ac, explicitly applying heavy quark 
symmetry to their trial wave functions. To date, there appear to be only two lattice studies of the semileptonic 
decays of heavy baryons. A first study of Ah and Sh semileptonic decays was made by Bowler et al. 1^. while 
Gottlieb and Tamhankar 0 have examined the decay of the Ah. Perez-Marcial and collaborators |l9| have 
studied the semileptonic decays of a number of charmed baryons, both in a non-relativistic quark model and 
in the MIT bag model. There have been light-front calculations |^, as well as ones using sum rules |^. 
Bethe-Salpeter formalisms [^, bag models [ 23 , and quark model calculations Large Nc arguments have 
also been applied to these form factors as well as perturbative QCD argu ments 0 . For the decays of 
a heavy baryon to a light one, work has been done using QCD sum rules |^, and there’s one quark model 
calculation apart from the work of Scora [^, to the best of our knowledge. 

The experimental status of heavy baryon semileptonic decays is somewhat rudimentary. The semileptonic 
decay rate for Ac ^ A has been measured by the CLEO and Argus collaborations , while the Delphi 

collaboration has only recently published an analysis of the exclusive semileptonic decay of the A^ |^. Prior 
to this, only the inclusive semileptonic branching fraction Af, ^ Acii^X had been reported in the PDG 0 . In 
their analysis of the Ac —> A semileptonic decay, the CLEO collaboration have assumed that the ground state 
A saturates the semileptonic decays of the Ac, and cite the absence of any final states of the form Aiv with 
additional decay products from the Ac to support their assumption |3l|. No experiments have yet reported 
results for the decay Ac —> n. 

The major difficulty in the baryon sector is that there is no source of heavy baryons as there is for mesons. 
Electron-positron colliders have produced billions of B mesons, utilizing the fact that the T(4s) is just above 
the BB threshold. In principle, a similar abundance can be duplicated among D mesons, by using the 'I'(3s). 
With baryons, production at such machines will be continuum production, as there are no (known) resonances 
to enhance the rate of production. Hadron colliders can provide larger yields, but they provide large yields 
of everything, and the heavy baryons will then have to be separated from everything else that is produced. 
However, the recent CLEO measurement suggests that some optimism regarding the future measurement of 
these decays might be warranted. In addition, there might be prospects for such studies at Jefferson Laboratory 
upgraded to 12 GeV or higher, or at E907 at FNAL. The advantage in these cases is that the target will be a 
baryon, unlike the continuum production of e“''e“ machines. 

In this paper we study the semileptonic decay of Ag baryons, the motivation for which is two-fold. One 
of our motivations is the importance of the CKM matrix elements Vub and 14b, and that baryon semileptonic 
decays can provide complementary extractions of these quantities, despite the difficulties mentioned above. In 
particular, a model such as ours, coupled with constraints provided by HQET, may lead to a more precise 
extraction of 14b than provided by meson decays. 

Our second motivation is to examine the predictions for these decays of a quark model developed very much 
in the spirit of the work by Capstick and Isgur , which builds on the work of Is gur and Karl 01 ■ Such a 
model has been applied, with some success, to the strong and electromagnetic couplings of baryons, 

and the semileptonic decays of baryons is a useful complementary extension of such a model. Indeed, a similar 
model, applied to the semileptonic decays of mesons 0, gave rise to HQET. We note that the thesis of Scora 
treats a number of baryon semileptonic decays in a framework very similar to that used in the treatment of 
mesons in 0. We use a similar framework, but we extend the model to examine the decays to excited baryons, 
whereas Scora examined only decays to ground state baryons. We also use a more sophisticated treatment 
of baryon structure. 

This manuscript is organized as follows: in Section H we discuss the hadronic matrix elements and decay 
rates. Section HI presents a brief outline of heavy quark effective theory as it relates to the decays that we 
discuss. In Section IV we describe the model we use to obtain the form factors, including some description of 
the Hamiltonian. Our analytic results are discussed in Section V, our numerical results are given in section VI, 
and Section VH presents our conclusions and outlook. A number of details of the calculation, including the 
explicit expressions for the form factors, are shown in a number of Appendices. 
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II. MATRIX ELEMENTS AND DECAY RATES 


A. Matrix Elements 


The transition matrix element for semileptonic decay of Ag (Ag ^ Aqiiyi) is 

T = - 75)M^,(Ag(p',s')|T^|Ag(p,s)), 


( 1 ) 


where GfIV^ = g‘^/{8M^) is the Fermi coupling constant, Mw is the intermediate vector boson mass, Vg^ 
is the CKM matrix element, and u77^(l ~ J 5 )u,^i is the lepton current. Since quarks are confined, the matrix 
element of the hadron current is described in terms of a number of form factors. We will build a model of the 
baryons we wish to study, and obtain approximations to the form factors that describe the hadronic matrix 
elements. For transitions between ground state {J^ = 1/2“'") baryons, the hadronic matrix elements of the 
vector and axial currents are 


(Ag(p ,s )|V/|Ag(p,s)) = u{p ,s) [Fi{q +F 2 {q ) 


Pm 


mAc 


F^iq^)- 




(Ag(p',s')|A;.|Ag(p,s)) = u{p', s') [ Gi{q )7^ + G2(g ) 


.,2\ Pm 


WAq 




mA„ 


u{p, s). 


75u(p, s). 


( 2 ) 

( 3 ) 


where the Fi and Gi’s are baryon form factors which depend on the square of the momentum transfer q = p — p' 
between the initial and the final baryons. Similarly, the matrix elements for decays to a daughter baryon with 
= 3/2“ are 


(A^^^(p',s')|V"^|Ag(p, s)) = u°‘{p',s') 
(A^^^(p',s')|A^|Ag(p, s)) = lt{p',s') 
The spinor u“(p', s') satisfies the conditions 


Po 


TOAq 

Pa 


TO^3/2 


ruA 


Fi-f^ + F2^+Fi 

^ mAQ 

Gi7^+G2^+G3 


Pu 


Q 


rriA, 


Q 


m^3/2 




p>“(p', s') = 0, U°‘{p', s')ja = 0, m“(p', s')/ = TO 3/2u“(/, s'). 
The corresponding matrix elements for decay to a baryon with = 5/2“'“ are 


u{p, s), 

75 m(p, s). (4) 

( 5 ) 


(Af (/,s')|y/Ag(p,s)) = ll“^(/,s') 
{^q^‘^ip\s')\A^\AQ{p,s)) = if‘^{p',s') 


TOAq 

Pa 


TOAq 


PP 

TOAq 

pp 

TOAq 


Pm 


Pm 


i^l7M + -P’2^^+i^3 


TOAc 


Tfl 


Gi7;x + G2^+G3 


WAq 


TO^5/2 


■ ^45/3m 

F Gigpfj, 


u{p, s), 
lbu{p, s), (6) 


where the spinor u°‘^{p'^ s') is symmetric in the indices a and /3, and satisfies 

Pau“^(p/s') = Ppu°'^{p,s')=0, 

{p', s')ja = m“^(/,s') 7^ = 0, 
U°‘^{p',s')f = TO^5/2 m“^(/, s'), 

U°'^ip',s')gaf3 = 0. 


( 7 ) 


Here we have only presented the form factor equations involving spinors having natural parity. The equations 
for unnatural parity spinors can be constructed in a similar manner by switching 75 from the equations defining 
the Gi to the equations defining the Fi. 
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B. Decay Rates 


The decay rate that arises from any of these matrix elements is 


(^n j - E^/), 

where A refers to the initial hadron. The leptonic tensor is 


The hadronic tensor is 


= 8[p>;p:^ • p., + ie^-'‘^f^pi^p.,0\. 




where Aq and refer to the initial and final baryons, respectively. The tensor Hf^i, must have the Lorentz 
structure 

= -agf,^+P++{p + p')f,{p + p')^+f3+_{p + p')f,{p-p')^ 

+ fi-+{p - p')tJ.iP + p')iy + P — {P-p')tJ.ip-p')iy 
+ ili^i^pa(p+p'Y(p-p'Y- 

The complete expression for the differential decay rate is 

(pY o G%m\ 

dYdYi = \yQ^\ P++Gp++ + P-+C0-+ + P+-C0+-+ P—C0..+^C^], ( 11 ) 


where 


= 


v-4-]^ 

^AqJ 


C'/3++ = 8 [x{2xm +y)- 2a:^ - y/Y\ - 


m'j ( mj 4to^ 


^Aq \^Aq ^Aq 


- — 8x + 3j/ 


C'/3_+ = C'/3+_ = 4(x - Xm) - y- 


jyL ( zl. 

2 \ y 2 

^Aq V 


El 

2 ’ 


= T‘2y 2xm,-4:X + y-\ - ^(2xm+y) ■ (12) 

L J 

In these expressions, x = EijmAQ, where Ei is the lepton energy, Xm = (^Aq ~ ”^a,)/(2w^^), and y = 
q^/m\^ = {p — p'Y The =F sign in Cry is determined by the charge of the lepton, with the upper 
(negative) sign corresponding to decays to £v(^. The lepton energy has the range 

A’ - Ai'l 




with K = 4[(2xm — yY — and y has the kinematic range niflrnjyy^ < 2 / < (^Aq ~ /'^\q- 

If the lepton mass is neglected, the terms in X ) P h P-- vanish, and the differential decay rate becomes 


= 


Q2 

P .r„ -^ + 2/3++[2x(2xm + 2 /) - 4a:^ - y] =F 7 y (2a;m - 4a; + 2 /) , (13) 
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where the lepton energy is now constrained by —it '/2 + {2xm + y)/4< x < Kl2 + (2xm + y)/4, and the lower 
limit on y is zero. In this case the differential rate depends only on a, /3++ and 7 . The explicit expressions for 
a, /3++ and 7 in terms of form factors for different final baryon spins are given in Appendix IDI 


III. HEAVY QUARK EFFECTIVE THEORY 

Heavy quark effective theory (HQET) has been a very useful tool in the study of electroweak decays of 
heavy hadrons. This effective theory has been applied to a number of processes, both inclusive and exclusive, to 
higher and higher order in the 1 /toq expansion, where mg is the mass of the heavy quark. In most applications, 
the aim has been to constrain the hadronic uncertainties in the extraction of CKM matrix elements such as Vub 
and Vcb- In this section, we take a different tack; we examine the predictions of HQET for decays of a heavy A 
into any of the allowed excited daughter baryons, whether this daughter baryon is heavy or light, with the aim 
of comparing these predictions with the form factors that we obtain in our model. 


A. Heavy to Heavy 


In a heavy excited baryon, the light quark system has some total angular momentum j, so that the total 
angular momentum of the baryon can be J = j ± 1/2. These two states are degenerate because of the heavy 
quark spin symmetry. It is useful to show explicitly the representation we use for these two degenerate baryons. 
In the notation of Falk ^3? we write {v') = {v') — with 

i^y') = (14) 

Here, uq{v) is the spinor of the heavy quark, and (r;') is a tensor that describes the spin-/ light quark 

system. This tensor is symmetric in all of its Lorentz indices, meaning that the {v') is also symmetric in 

all its Lorentz indices. Both (I'O satisfy the conditions 


= 0, (15) 

where /i^ and ni indicate any pair of the indices gi... gj. The state with J = j -|- 1/2 also satisfies 


7/x.W^-+1/2 


• 1^3 


= 0 . 


(16) 


Further details of the structure and properties of these tensors are given in Falk’s article [4(T| . 

At this point, it is useful to discuss the parity of the states, which is determined by the parity of the light 
component. A spin-j light quark component with parity (—1)^ is said to have ‘natural’ parity, unnatural parity 
otherwise. The natural-parity light quark systems therefore have = {2n)'^ or = (2n + 1)“, with n a 
positive integer or zero. The natural-parity light quark systems are represented by tensors, while those with 
unnatural parity are represented by pseudo-tensors. Since the parity of the baryon is that of the light quark 
system, we may refer to the baryons as being tensors or pseudo-tensors, with the understanding that this really 
refers to the light-quark component of the baryon. It is thus convenient to divide the decays we discuss into 
two classes, those in which the daughter baryons are tensors, and those in which they are pseudo-tensors. We 
begin with the discussion of the tensor decays. 

In general, we are interested in the matrix element 


A =< A*(w',j)|cr6|A{,(u) >, 


(17) 


where c and b are the heavy quark fields, and T is an arbitrary combination of Dirac matrices. With the use of 
HQET, we may write this as 

< A:(u',j)|cr6|Ab(u) (18) 

to leading order. In writing this form, we are omitting multiplicative QCD corrections of order unity that arise 
from matching of the effective theory to full QCD at different mass scales. Here, is the most general 
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tensor that we can construct, given the kinematic variables at our disposal. Clearly, may not contain 

any factors of or , and therefore takes the form 

='q^^\v-v')vf,^...v^^. (19) 

Thus, a single form factor, ■ v') is needed to this order, regardless of the spin of the final baryon. In 

addition, spin symmetry allows us to relate the form factors for T = 7 ^ to those for T = 7 ^ 75 . 

The case of = 1/2+, j = 0 requires a special comment. These states may be thought of as radial 
excitations of the ground state A+. Because of the heavy quark symmetry, and the orthogonality of these states 
with respect to the ground state, we must have 

< = 0+)(„)|cr&|Ab(z;) >= {vv - • v)uiv)Tu{v), (20) 

where the subscripts (n) denote the nth radial excitation. That is, these amplitudes must vanish as v' v. 
This result has been pointed out by Isgur, Wise and Youssefmir Note, too, that all of the other amplitudes 
(j 7 ^ 0 ) vanish trivially at the non-recoil point. 

For the pseudo-tensor decays, we write exactly the same form, but must now be a pseudo-tensor 

object, and must therefore be constructed by using the £ tensor. Inspection shows that no such pseudo-tensor 
can be constructed, given that we have only two kinematic variables at our disposal, namely v and v', and that 
the spinor-tensor used to describe the daughter baryon is symmetric in its indices. Thus, decay amplitudes for 
transitions to pseudo-tensor daughter baryons vanish at leading order in HQET. 

Applying these results to the specific case of = 1“, we find, for = l/2“, 

Fi = F2=G2 = G3 = F3 = 0, Gi = (w). (21) 


For 3/2-, 


F 2 = F 3 = G 2 = G 3 = F 4 = G 4 = 0, Fi = Gi = ( 22 ) 

In these two sets of equations is a universal function of the Isgur-Wise type, and w = v ■ v'. 

For = 2+, we find for = 3/2+, 

= G 3 = F 4 = G 4 = 0, F2 = G2 = Gi = (23) 

and for 5/2+ 


F 2 = F 3 = F 4 = G 4 = G 2 = G 3 = 0, Fi = Gi = (24) 

As with the previous example, the function is an Isgur-Wise form factor common to both decays. 

For the elastic decays, as well as for decays to the 1/2“, 3/2“ doublet, the matrix elements have been 
evaluated at order 1 /rric and 1 /nib in the heavy quark expansion . When we present our results for the form 
factors, we will compare our expressions with the predictions of HQET. 


B. Heavy to Light 

For the heavy to light transitions, we may no longer describe the daughter baryons in terms of the spin 
structure of the light quark system that helps to make up the baryon. Instead, we are forced to use the total 
angular momentum of the baryon concerned, as well as its parity. As before, we may represent one of these 
baryons, denoted A*, by a generalized Rarita-Schwinger field {p), where the auxiliary conditions now are 

^yMi. ./^r.(p) ^ mA*M+^' +"(p), 7 ;,iM'"^- +"(p) = 0, 

«'"'■■+"(?) = 0, u^'" '^"(p) = 0, 


( 25 ) 







and a baryon with angular momentum and parity is represented by a spinor-tensor with n = J — 1/2 indices. 
As was the case with the heavy to heavy transitions, we need to divide the possible transitions into two classes, 
which we call tensor and pseudo-tensor, with the obvious meaning. 

As before, we begin with the transitions to tensor states. Here, we say a state of total angular momentum J 
is a tensor if its parity is (—and is a pseudo-tensor otherwise. The matrix element of interest is 

< A*(p)jp|src|A+(u) >= (26) 

where is the most general tensor that one can construct, and n = J — 1/2. As with the heavy to heavy 

transitions, we may not use any factors of 7 ^^, or in constructing which must therefore have 

the form 


.. .Vfj.^An- (27) 

Here, An is the most general Lorentz scalar that we can build. On inspection, we find that 

An=£.t^ (28) 

so that each of these transitions is described by two form factors, at leading order in HQET. 

For the transitions into pseudo-tensor daughter baryons, we write exactly the same form as in Eq. 
but now must be a pseudo-tensor. This may involve the use of the e tensor, but since u^i -^"(p) 

is symmetric in its indices, at most one of these indices may be contracted with the indices of the e tensor. 
With some patience, and the use of a few well chosen identities, one can show that any pseudo-tensor term 
constructed with the e tensor may always be reduced to an ordinary tensor multiplying a 75 matrix. We will 
therefore leave out much of the tedium, and simply write for these transitions 

-h 75 , (29) 

where the and (/i are functions of the kinematic variable v ■ p'. Thus, any of the heavy to light transitions is 
described by a pair of form factors, to this order in HQET. Note that for both sets of heavy to light transitions, 
we may use the spin symmetry of HQET to relate the two form factors necessary for T = 7 ^ to those for 

r = 7 m75- 

For 1/2+, we find 


F 3 = G 3 = 0 , F2 = G2 = 2 d°\ El = Gi = 

while for 1/2“, the form factors are 

F3 = G3 = 0, F2=G2 = -2d°\ El = - (cr + , Gi = - ■ 


For 3 / 2 “, 

E3 = G3 = E4 = G4 = 0, E2 = G2 = 2d'\ El = Gi = 

For 3 / 2 +, 

E3 = G3 = E4 = G4 = 0, E2 = G2 = -2cd\ Fi = - (d'^ + cd^), Gi = - (d'^ - d'’) ■ 

For 5 / 2 +, 

E3 = G3 = E4 = G4 = 0, E2 = G2 = 2d"\ El=d"^ - d"\ Gi = d'^+d'^- 


(30) 

(31) 

(32) 

(33) 

(34) 


Note that, in principle, the form factors for the decays to 1/2 have no relationship with those for decays to 
3/2“, in this limit. 
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IV. THE MODEL 


A. Wave Function Components 


Our calculation follows the spirit of the work by ISGW 0. In our model, a baryon state has the form 


I^q(P. s)) = j '^^PpC^^?’AC'^^AQkl(Pl:Sl)'?2(P2,S2)93(P3,S3)), 

where Pp = ■;^(Pi — P2), Pa = "^(Pi +P2 ~ 2p3) are the Jacobi momenta, is the antisymmetric color wave 
function and ^Aq = 4’Aq'^AqXAq is a symmetric combination of flavor, momentum and spin wave functions. 
For Aq the flavor wave function we use is 


'Paq 


V2 


{ud — du)Q, 


which is antisymmetric in quarks 1 and 2. The momentum-spin portion of the wave function must therefore 
be an tisy mmetric in quarks 1 and 2. For states like the neutron and proton, we use the ’’uds' basis used in 
Refs. [j.Jl l34| . In that basis, the wave function of the proton is simply uud, while that for the neutron is ddu. 
This flavor wave function provides some simplification in dealing with matrix elements of the Hamiltonian. 
However, the treatment of current matrix elements, such as those that describe semileptonic decays, will require 
some extra care, as will be explained later. 

The total spin of the three spin-1/2 quarks can be either 3/2 or 1/2. The spin wave functions for the maximally 
stretched state in each case are 


xf/2(+3/2) = ITTT), 

x//,(+i/ 2 ) = ^(itiT)-UTT)), 

x^/2(+i/ 2) = -^(|tiT) + UTT)-2|m)), 


where S labels the state as totally symmetric, while A(p) denotes the mixed symmetric states that are symmetric 
(anti-symmetric) under the exchange of quarks 1 and 2. The momentum wave function for total L = ip + £x 
is constructed from a Clebsch-Gordan sum of the wave functions of the two Jacobi coordinates Pp and pa, and 
takes the form 


'fpLMnpl.n^ixiPp^Px) = pTU, i^M - -m{P\) ■ 

m 

The momentum and spin wave functions are then coupled to give symmetric wave functions corresponding to 
total spin J and parity (— 1 )(V+^a)^ 

^JM = SM — Mi,)'lpLMLnpEpnxixXs{M — M^). 

Ml 


The full wave function for a state A is built from a linear superposition of such components as 

'^A.J-M=^AY.’nt'^^JM- (35) 

i 

Here ^a is the flavor wave function of the state A, and the rjf are coefficients that are determined by diagonalizing 
a Hamiltonian in the basis of the 'I'jm- For this calculation, we limit the expansion in the last equation to 
components that satisfy N < 2, where N = 2(np + n\) + £p + £\. Consistent with this is the fact that the states 
we discuss all correspond to V < 2. With this limitation, the wave function for a Ag with = 1/2+ takes 
the form 
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'I'Aq,i/ 2 +m — ?^>Aq I ? 7 f‘^V'oooooo(Pp, Pa) + '0001000(Pp,P a) + ?7^®'0ooooio(Pp,P a) Xi/ 2 (-^) 


0000101 (Pp, Pa)XV 2(^) + % ® 01Mi,OlOl(Pp, Pa)X3/2(^ - Ml) 


+ % 


01 MiOlOl(Pp,PA)XV 2 (^ - ^a) 


J 1/2,M 


Aq 

^7 


1/2,M 
S 


(36) 


02 MpOlOl(Pp,PA)xi/ 2 (^ - ^a) 


1/2,M 


where [0LMpnp^pn^^;, (Pp, PA)xs(Af - Ml)] is a shorthand notation that denotes the Clebsch-Gordan sum 
J2ml^'^M\LMl,SM - ML)iliLMLn^e^n^e^{Pp,P\)xs{M - Ml). When we diagonalize the Hamiltonian, this 
expansion will provide the wave functions for seven states with = 1 / 2 +, the lowest of which will be taken 
to be the ground state of the system. 

A simplihed version of the model would truncate this expansion after the first component, giving 


'I'Aq, 1/2+M — 0A0OOOOOo(Pp, Pa)x^/ 2 (-^)) 
while the first radial excitation of interest in this model would be 


'^Aq,i/2+m — 0a0ooooio(Pp,Pa)x?/2(-^)- 
There exists a second radial excitation which, in the truncated basis would be 


'I^Aq,1/2+M — 0A0001000 (Pp,Pa)x)’/2(-^)- 


The latter state has its radial excitation in the p coordinate, which means that it has a very small overlap 
with the ground state in the spectator model that we use. For some states, this truncation provides a very 
good approximation to the wave function, but there are important conhguration mixing effects for a number 
of states. In the spectator assumption that we use, not all of these states have an overlap with the initial 
ground-state Aq. The possible states which can be connected to the ground state are the states with = 
1/2+, 1/2'^, l/2“, 3/2“, 3/2+, 5/2+, where 1/2+ and 1/2+ denote the ground state and the hrst (radially) excited 
state. 

It is useful for us to list the single-component representations of these states. The states with J+’ = 1/2+ 
have already been given. For the remaining states, we have 


'^Aq,1/2-M 

^Aq,3/2-M 

'^'Aq,3/2+M 

'^'Aq,5/2+M 


0A 

0A 

0A 

0A 


01 MiOOOl(Pp,PA)Xi/ 2 (^ 
01 MiOOOl(Pp,PA)x)’/ 2 (-^ 
02MiOOO2(Pp, P\)Xi/2iM 

02Mi 0002 (Pp, PA)x^/2 (-^ 


Ml) 

Ml) 

Ml) 

Ml) 


1/2,M ’ 
3/2,M ’ 
3/2,M ’ 
5/2,M 


(37) 


From these representations, the multiplet structure expected in the heavy quark limit is easily identified, with 
the 1/2“ and 3/2“ states forming a multiplet, and the 3/2+ and 5/2+ states forming another. Both of the 
1 / 2 + states we consider are singlets. 


1. Expansion Bases 

A common choice for constructing baryon wave function is the harmonic oscillator basis. One advantage of 
using this basis is that it facilitates calculation of the required matrix elements. However, it leads to form factors 
that fall off too rapidly at large values of momentum transfer. We therefore also use the so-called Sturmian 
basis |42| . In this basis, form factors have multipole dependence on , which is what is expected experimentally. 
The full wave functions in momentum space are 
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^nLm{V>) = 


2n\ 


(n + L + i)! 

in the harmonic oscillator basis, and 


(i)^(-l)"'—^-3-6 (p^/a^)3^Lm(p) 

a^+i 


(38) 


,.St / X 2[n!(n + 2L + 2)!]2 , 1 

V^nLm(p) / I r I l\i 


(n + L+i)! + 


L +2 




^L+|,L+i^ 


2’^+2 J p^-13- 


+ /32 


3^Lm(p) (39) 


in the Sturmian basis. The L'^{x) are generalized Laguerre polynomials and the {y) are Jacobi polynomials, 
with p = |p|. The corresponding wave functions in coordinate space are 


V’nLm(l') = 

in the harmonic oscillator basis, and 
V’fLm(l’) = 

in the Sturmian basis. 


2 n! 


(n + L + i)! 


a^+2e “s" Ln^‘^{a'^r‘^)yLrn{'r) 


_(n + 2 L + 2 )! 


1 

" (2/3)^+ie-^’'Lf+2(2/3r):yL™(r) 


2 . Hamiltonian 


We use a non-relativistic quark model similar to that of Isgur and Karl I^H^, with some of the modifications 
suggested by Capstick and Isgur The Isgur-Karl model evolved from the pioneering work of others; an 

extensive list of references to the origins of the model can be found in Ref. . 

The phenomenological Hamiltonian we use takes the form 

^ = E^*+E(o+<p). (40) 

i i<j 

where Ki is the kinetic part of the Hamiltonian. For this, we use two forms, the usual non-relativistic form 
given by 

= (”■• + i-) ■ 

and a semirelativistic form given by 


Ki = ypf + ^l 

The spin independent confining potential is a simplified version of that used by Capstick and Isgur 


TAb ^ 1 2q!Cou1 

*^conf ~'-^999 ■*" 2 3 r ■ ' 


with rij = \ri — rj|. Here f?hyp is the hyperfine interaction, assumed to have the form 

- s,; • S, 


HZ^ = <! ^ 


irriim^ \ 3 




3(Si • Ty )(Sj • Tij) 


(42) 
I, with 

(43) 

(44) 


The first term is a contact term, while the second is a tensor term. The spin-orbit interaction is neglected. We 
note here that acoui, cthyp, b, Cqqq, and rrii are not fundamental, but are phenomenological parameters obtained 
from a fit to the spectrum of baryon states. 
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B. Obtaining the Form Factors 


1- Aq Aq 


Here, we illustrate the procedure we follow to obtain the form factors, using the decay of the Aq to the ground 
state Aq as an example. We begin with the vector current matrix element from Eq. with the assumption 
that the parent Aq is at rest and the daughter Aq has three momentum p. The left-hand side of Eq. lO is 
evaluated using the quark model, after the operator has been reduced to its Pauli (non-relativistic) 

form. Specific values for the index /r are chosen, as well as specific values of s and s'. By making three sets 
of such choices, three equations for the Fi in terms of the quark-model matrix elements of three operators are 
obtained. This system of equations is then solved to obtain the expressions for the form factors. In the specific 
case at hand, choosing s = s' = +1/2 and /r = 0, for instance, leads to 


where 


(A,(p,+)|g7oQ|AQ(0,+)) = J d^p'^d^p'^d^ppd^pxC^*C^'+f^{+) 


X {qW2Q\<ihoQ\qiq2Q)'^AQi+) 

= Fi + +2 + Fs, 


(9Wkk^7oQ|gig2Q) = {qW2\QiQ2){q\qhoQ\Q)- 


(45) 


(46) 


The matrix element (qi 92 I'll 92 ) gives ^-functions in spin, momentum and flavor in the spectator approximation, 
while the operator 970 Q = 1 + 0 { ^ ^ ). Using the ^-functions, the integral is simplified to 

\ 7nqTri(^ j 


A9(p,+) 


1 + 0 


1 


mqlTiQ J 


AQ(0,-k))= / d^ppd^px'i/liv'p.v'x) 


1 + 0 


1 


TTlqmQ 


V'Ao(Pp,Pa), (47) 


with Pp = Pp, Pa = Pa — 2^/372 mcrp/wA,) where rria is the mass of the light quark. This leaves the momentum 
integration, which is performed by using both bases for the momentum wave function shown earlier. The analytic 
results for the form factors for Aq decaying into various A^ final states are given in Appendix El For decays to 
excited states, the calculation of the form factors is a little more involved, but the basic idea is as outlined here. 


g. Aq 

For decays in which the daughter baryon is a nucleon, the procedure is much the same as outlined in the 
previous subsection, with one modification. To illustrate, let us take the specific example of A^ —> p. The flavor 
wave functions of these two states have been chosen to be 

(pAt = —^{ud— du)b, (j)p = uud. (48) 

v 2 

For the transition to occur, the third quark in the parent baryon, the b quark, undergoes the transition b u, 
leaving a final state that is '^{ud — du)u. This has no overlap with the flavor wave function that we use for 
the proton. We must now permute the third quark with the first and second quarks, giving 

{IZ}—j=(ud — du)u =—j={udu — uud) ^ {23}—^{ud — du)u=—^{uud — duu), (49) 

v 2 v 2 v 2 v 2 

both of which now have some overlap with the proton flavor wave function we use. This requires that the sum 
of matrix elements 


(A^(p, -l-)|{13}Oi|AQ(0, -I-)) -I- (A^(p, -|-)|{23}Oi|AQ(0, -I-)) 
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be evaluated, where we apply the permutation to the wave function of the daughter nucleon. The permutation 
operators also transform the spin and momentum wave function of the nucleon. The transformed spin wave 
functions are 


{13}x^(s) = - ^X^(s), {23}x^(s) = ^X^(s) - 


(50) 


After carrying out the transformation on the nucleon wave function, and using the fact that the ground state 
momentum space wave function is totally symmetric, we find 


(p(p,s)|O,|Aq(0,s')) = (-vW4) / d^Ppd^p\i)*p{v'p,v'\)A""'{Oi)i}AQ{Vp,V\), 


(51) 


where A®'* {Oi) is the Pauli reduction of the operator Oi. The integrations required for the Af, to proton form 
factors are the same as those in Eq. (in in the previous subsection, and so the form factors are the same up 
to a multiplicative factor. For excited states, however, the procedure is slightly more involved, and is easily 
illustrated by examining the decays to the radially excited nucleon. 

Assuming single components, the wave function of the radially excited state is 


^Ar,l/2+M = '('iV'0OOOOlo(Pp,PA)Xl/2(^)- 

The {13} transformation, acting on the spin-space part of this wave function, produces 


{13}'I'ooooio(Pp, Pa)Xi/ 2 (-^) 


'l'ooooio(Pp, Pa) 


2 






= v^V'ooiooo(Pp,Pa)xv2(^)+3'!/'ooiooo(Pp,Pa)xv 2(^) 
-h y3V'ooooio(Pp,PA)x)'/2W + V'ooooio(Pp,Pa)Xi/2W 
-f yT8V'00010l(Pp',PA)x)'/2W + v^'0oooioi(Pp,Pa)xv2(^) : 


(52) 


(53) 


with a similar expression for the {23} transformation. Here Pp' =-^(pa — P 2 ), Pa' = -^(Ps + P 2 — 2pi) are the 
Jacobi coordinates in the transformed basis. Of these components, only the first, third and fifth have spin wave 
functions that overlap with the decaying Aq, while only the first and third have non-zero spatial overlaps. The 
integrals that arise from the first component are simply a numerical factor (\/27/8) times those that arise in 
the Aq — > Kq matrix elements, for the radially excited Ag. The integrals that arise from the third term are also 
a numerical factor (a/3/8) times the Aq ^ A, ground-state integrals, multiplied by a factor that arises from 
the spectator overlap. In this case, this overlap is expected to be small, since the spectators are in a radially 
excited state in the daughter baryon, but in their ground state in the parent. 

The above procedure is relatively straightforward to implement in the harmonic oscillator basis, largely due 
to the fact that the Moshinsky rotations have been treated by a number of authors, and are also fairly simple 
to calculate. In particular, the fact that the ‘permuted’ wave function can be written in terms of a finite set of 
transformed wave function components is another feature that makes the harmonic oscillator basis attractive 
for calculations like these. In the Sturmian basis, however, the permutation of particles requires an infinite 
sum of transformed wave functions. This sum could be truncated at some point in a calculation such as this. 
However, at this point we do not examine decays to daughter nucleons in the Sturmian basis. 


V. ANALYTIC RESULTS AND COMPARISON WITH HQET 

The analytic expressions that we obtain for the form factors are shown in AppendixEl for both the Sturmian 
and harmonic oscillator bases. The results shown there are valid when the wave function for a particular state 
is written as a single component, in either expansion basis. 

As mentioned earlier, one of the advantages of the Sturmian basis is that it leads to form factors that behave 
like multipoles in the kinematic variable, and this is seen in the forms that we display. At this point, it is 
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instructive to compare, as far as possible, these analytic forms with the predictions of HQET. While HQET 
does not give the explicit forms of the form factors, a number of relationships among the form factors are 
expected, and any model should reproduce these relationships. In what follows, we restrict our comparison to 
the predictions that are valid at the non-recoil point, as we have ignored any kinematic dependence beyond the 
Gaussian or multipole factors shown in Appendix [0 In addition, we focus mainly on the predictions for heavy 
to heavy transitions. 


A. Natural Parity Daughter Baryons 


We begin by discussing the form factors for decays to daughter baryons of natural parity. In this work, this 
means daughter baryons with = 1/2+ (both ground state and first excited state), = 1/2“ and 3/2“ 
(which constitute a degenerate doublet when the daughter baryons are also treated as heavy) and = 3/2+ 
and 5/2+ (also a doublet). In our discussion of these results, we implicitly assume that the wave functions for 
the states are dominated by a single component of the wave function expansions that we use. These single¬ 
component wave functions have been described in section HV Al 

For elastic decays, predictions have been made at least to order 1/m^ and 1 /toq. However, we will restrict 
our discussion to the predictions valid to order Ijrriq and l/mg. To this order, using the results of Falk and 
Neubert ^5, the relationships among form factors are 


F2 

Fi 


(-T2 — - rs — -Da, 


Too 


Too 


Gi — F 2 



(54) 


Our expressions for the form factors satisfy these relationships, in both bases, to the appropriate order. In fact, 
the analytic forms obtained exactly match the structure predicted by HQET . 

For the (l/2“, 3/2“) doublet, there are 14 form factors in general, which Leibovitch and Stewart M write in 
terms of a number of universal functions and constants, valid at order Ijmq and XjmQ. Using their expressions, 
and writing form factors for the 1/2“ state as primed quantities, the relationships expected are 


F[ 


^3 

Fz 

Gi 


U' + G' + 2G; 


2V3n 


(3toq - mq) Fi, 


3G' + ^ (G 3 - 2 F 4 ), 

—G3, G2 — F2, F\ — G2 = G3 — F2, 

- 3 F 4 + 2 V 3 G', E'-G'=-^G 3 , 


ViFi { 1 -k 


mg 




(55) 


where terms that vanish at the non-recoil point have been ignored. Our results for these states also satisfy all 
eight of the relationships shown above, in both bases. Thus, there is a very good correspondence between the 
predictions of HQET and those of the quark model that we use, and this correspondence is independent of the 
wave function basis chosen. 

For the (3/2+,5/2+) doublet, the available predictions are at leading order, shown in Eqs. 123) and 
These are also satisfied by our analytic expressions for the form factors, in both bases. 

For the excited state with J+ = 1/2)^, the predictions of HQET are that the form factors should vanish at 
the non-recoil point, by reason of the orthogonality of the wave functions. In the treatments in the literature, 
this is achieved by assuming that the form factors have an explicit factor that vanishes as w —> 1. In the 
expressions that we have obtained for the leading order form factors, this orthogonality arises explicitly from 
the size parameters of the wave functions. 

It is instructive to examine the expression for Fi for this decay, in the limit when the Hamiltonian is that of 
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a harmonic oscillator. The expression for Fi is 


Fi=Ih 


1 




{al-aD- 


nio 


iaXx' 


a% 


mq niq 


where 


/ 3/2 3/2' 



a' 


exp - 


AA' 


3m^ 

2wL a^A' 


(56) 


(57) 


In the above expressions, aA(crA') is the size parameter of the initial (final) wave function associated with the 
Jacobi coordinate A, and a\x, = (a a + c«a')/ 2- If the Hamiltonian is taken to be a harmonic oscillator of the 
form 


K 


(^|ri - r2|^ + |ri - + |r2 - = ZK + A^) 


h" = y ( Fi - 1-2 


(58) 


where is the position of the z-th quark and p = (ri — and A = (ri + r 2 — 2r^)/x/Q are the Jacobi 

coordinates, then 


/ iKmaTUQ / “^KmaTriq \ 

VTOQ + 2mcr/ ’ \mqF2ma) 


(59) 


With these forms, the term in Fi proportional to ma vanishes identically, while the term in {a\ — a\,) becomes 
proportional to \lmq — IjraQ, and so vanishes in the heavy quark limit. The terms in which we do not 
include here, will be those that contribute, despite the orthogonality of the wave functions, as expected. Note 
that even though the terms will appear with explicit factors of l/rn^, p will range from small values (of order 
Aqcd)j to a maximum of (w^^ — m\^)/{2mAQ). Such terms are therefore not necessarily negligible. However, 
in the non-relativistic model that we use for the form factors, we have neglected such terms. 


B. Unnatural Parity Daughter Baryons 


For the decays to baryons with unnatural parity, HQET predicts that the form factors should vanish at 
leading order. In the present model, we first have to identify such states, which we do in the heavy quark limit, 
using the single-component wave functions. The wave functions of interest are 


^Aq,1/2+M 

'I'Aq,3/2+M 

'I'Aq,3/2-M 

'I'Aq,5/2-M 


(t>A i^oooioiiPp,Px)Xi/2iM - Ml) 
(t>A [V'00010l(Pp,PA)xf/2(Af- AIl) 


J 1/2,M 


J 3/2,M 




V'iMz,oioo(Pp, P\)xt/ 2 {M - Ml) 


(t)A il^iMLOiooiPp,Px)xi/2iM - Ml) 


3/2,M ' 
5/2,M ' 


(60) 


In the spectator assumption that we use, none of these states have any overlap with the ground state parent 
Aq. In fact, there is a ‘two-fold’ orthogonality at play. The spin wave function of the two spectator quarks 
is orthogonal to the corresponding wave function in the parent baryon. The spatial wave functions of these 
two quarks are also orthogonal in parent and daughter. Thus, decays to these states will only occur through 
configuration mixing in the wave function, induced by various terms in the Hamiltonian. 

In the model that we use, configuration mixing in the spin wave functions arises from hyperfine terms involving 
the heavy quark, which means that such mixing will be small. Thus we expect that decays to such states should 
be significantly suppressed. Interestingly, the suppression of the decays to these unnatural parity doublets 
persists as the mass of the heavy quark in the daughter baryon is decreased, as such configuration mixing 
remains small. In this case, even though the definition of unnatural parity is different for light states, there are 
still a number of decays (in Ac ^ A, for instance) that are predicted to be significantly suppressed. We will 
comment on this later, when we examine the numerical results of our model. 
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VI. NUMERICAL RESULTS 
A. Model Parameters, Mass Spectra and Wave Functions 

In Section llV A 21 we introduced the two Hamiltonians we diagonalize to obtain the baryon spectrum. The 
two Hamiltonians differed only in the form chosen for the kinetic portion, one of which was nonrelativistic (NR), 
while the other was semirelativistic (SR). In addition, we use two different expansion bases to obtain the wave 
functions: the harmonic oscillator (HO) basis, and the Sturmian (ST) basis. In the following, the four spectra 
we obtain will be denoted HONR, HOSR, STNR and STSR, in what should be an obvious notation. 

There are eight free parameters to be obtained for each spectrum: four quark masses (m„ = m^, mg, nic and 
mb), and 4 parameters of the potential (ohyp, acoui, and Cqqq). We have investigated the effects of a tensor 
interaction in the two harmonic oscillator models, and found the effects to be small. In the results we present, 
the tensor interaction has therefore been ignored. The eight parameters are determined from a ‘variational 
diagonalization’ of the Hamiltonian. The variational parameters are the size parameters Op and a\ of Eq. 
or f3p and /3 a of Eq. (PI- This variational diagonalization is accompanied by a fit to the known spectrum. In 
this fit, the eight parameters mentioned before are varied. The values we obtain for the Hamiltonian parameters 
are shown in Table ^ while some of the wave function size parameters are shown in Table im 


TABLE I: Hamiltonian parameters obtained from the four different fits. In the first column, HO refers to the harmonic 
oscillator basis, while ST refers to the Sturmian basis. In the same column, NR indicates a non-relativistic Hamiltonian, 
while SR indicates a semirelativistic one. The form of the Hamiltonian is described in Section fiV A 21 


model 

(GeV) 

rus (GeV) 

me (GeV) 

mt (GeV) b (GeV^) 

OCoul 

Ohyp 

Cqqq (GeV) 

HONR 

0.40 

0.65 

1.89 

5.28 

0.14 

0.45 

0.81 

-1.20 

HOSR 

0.38 

0.59 

1.83 

5.17 

0.17 

0.09 

0.26 

-1.45 

STNR 

0.40 

0.64 

1.87 

5.28 

0.13 

0.35 

0.31 

-1.22 

STSR 

0.34 

0.57 

1.78 

5.22 

0.15 

0.19 

0.11 

-1.23 


We note that the value of 6, the slope of the linear potential, tends to be smaller than in most published studies 
of the baryon spectrum. The same is true for the strength of the hyperfine interaction, ahyp- In the case of the 
latter, the small strength arises because the hyperfine interaction is treated as a contact interaction, and this can 
lead to very strong attractive forces between the quarks. One result of this is that, for sufficiently large values 
of Ohyp, the masses of the lightest baryon states can become negative. The small value of this parameter that 
results from our fits is therefore driven largely by the need for positive baryon masses. One direct consequence 
is that hyperfine splittings are not well reproduced in all but the HONR model, with the A — N mass splitting 
being about one third of its experimental value. 


In general, we allow the values of Op to be different from a\. The exceptions occur in cases when the three 
quarks are identical, as they are in the nucleon. In that case, the variational diagonalization automatically 
selects ap = a\. In Table m we show only some values of the size parameters. The other size parameters, for 
the states that are significant for this work, are related to those presented. For instance, for the 1/2)'' states, 
the size parameters are the same as for the 1/2''' states. Furthermore, since we do not include a spin-orbit 
interaction in our Hamiltonian, the size parameters for the 1/2“ and 3/2“ states are identical. We do not show 
the size parameters for the Aq states with Q = b, c, oi s and = 3/2''' or 5/2''', mainly because we find that 
semileptonic decays to these states are very small. 


1. Mass Spectra 


Portions of the four mass spectra we obtain are shown in Table nni In this table, the first two columns 
identify the state and its experimental mass, while the next four columns show the masses that result from the 
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TABLE II: Wave function size parameters, ttp and a\, for states of different in different models. All values are in 
GeV. For the Sturmian basis, the size parameters have been denoted /3 in the text. 



model 

A(. A, AN 

{a\, Op) (oa, Op) {a\, oip) (oa, Op) 

1/2+ 

HONR 

(0.59, 0.61) (0.55, 0.58) (0.49, 0.53) 0.48 

1/2+ 

HOSR 

(0.68, 0.68) (0.60, 0.61) (0.52, 0.57) 0.54 

1/2+ 

STNR 

(0.44, 0.66) (0.41, 0.69) (0.35, 0.75) 

1/2+ 

STSR 

(0.46, 0.64) (0.43, 0.67) (0.38, 0.72) 

1/2- 

HONR 

(0.47, 0.49) (0.40, 0.47) 0.37 

1/2- 

HOSR 

(0.55, 0.59) (0.48, 0.54) 0.46 

1/2- 

STNR 

(0.60, 0.50) (0.55, 0.54) 

1/2- 

STSR 

(0.61, 0.49) (0.58, 0.51) 

3/2+ 

HONR 

0.35 

3/2+ 

HOSR 

0.44 

5/2+ 

HONR 

0.35 

5/2+ 

HOSR 

0.46 


models that we use. The small hyperfine interaction that we alluded to in the previous subsection has resulted 
in ground state nucleons that are too heavy, in all models. In addition, the ground state A (not shown in the 
table) is too light in all models. Similar patterns emerge when the various Ag and Sq (not shown) states are 
compared. The size of this interaction also results in ‘radial’ excitations that are too heavy, even heavier than 
usually result in models like these. 

We note, too, that the different models give very similar results for many of the states such as the A(l/2+), 
N{l/2~), A(l/2+) and Af,(l/2+), for instance, but for some states such as A(l/2)'")(1440), there are striking 
differences in the masses obtained. 


TABLE III: Baryon masses in GeV fitted in different quark models. The first two columns identify the state and its 
experimental mass, while the next four columns show the masses that result from the models that we use. 


State 

Experimental Mass 

HONR HOSR STNR STSR 

A(l/2+) 

0.94 

1.00 

1.08 

1.08 

1.08 

iV(l/2+) 

1.44 

1.76 

1.60 

1.81 

1.70 

iV(l/2-) 

1.54 

1.45 

1.44 

1.50 

1.47 

A(3/2-) 

1.52 

1.45 

1.44 

1.50 

1.47 

A(3/2+) 

1.72 

1.72 

1.69 

1.78 

1.77 

A(5/2+) 

1.68 

1.72 

1.69 

1.78 

1.77 

A(l/2+) 

1.12 

1.23 

1.23 

1.12 

1.10 

A(l/2+) 

1.60 

1.73 

1.81 

1.61 

1.55 

A(l/2-) 

1.41 

1.54 

1.62 

1.50 

1.56 

A(3/2-) 

1.52 

1.54 

1.62 

1.50 

1.56 

A(3/2+) 

1.89 

1.81 

1.81 

1.77 

1.87 

A(5/2+) 

1.82 

1.82 

1.81 

1.77 

1.87 

Ac(l/2+) 

2.28 

2.35 

2.32 

2.26 

2.22 

Ac(l/2-) 

2.59 

2.61 

2.70 

2.61 

2.68 

Ac(3/2-) 

2.63 

2.61 

2.70 

2.61 

2.68 

A6(1/2+) 

5.62 

5.62 

5.62 

5.62 

5.62 
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Z. Wave Functions 

For many of the states that we treat, the wave functions that result are, to a very good approximation, the 
single component wave functions shown in Section flV Al This turns out to be a particularly good approximation 
for the orbitally excited states such as the 1/2“ and 3/2“ states, for all but the nucleon states. For the A(l/2“) 
and A(3/2“), for instance, the dominant component has a coefficient [the r]i of Eq. IIA2II ] of at least 0.985 in all 
of the models. We treat such states as being single component states, and this will introduce errors of about 
a few percent (typically less than three percent for the particular states mentioned, usually much less for the 
states containing a c or 6 quark). 


TABLE IV: Mixing coefficients {rji) of the two lowest lying 1/2"*' states in different flavor sectors. The r/i are defined in 
Eq. m of Appendix m 


Baryon states 

hi 

HONR 

V2 

m 

hi 

HOSR 

h2 

h3 

hi 

STNR 

h2 

h3 

hi 

STSR 

h2 

h3 

lV(l/2+) 

0.979 

-0.150 

0.034 

0.989 

-0.110 

0.028 

- 

- 

- 

- 

- 

- 

lV(l/2+) 

0.022 

0.522 

0.825 

-0.026 

0.579 

0.800 

- 

- 

- 

- 

- 

- 

A(l/2+) 

0.994 

0.005 

-0.069 

0.998 

0.003 

-0.035 

0.900 

0.208 

0.382 

0.875 

0.313 

0.368 

A(l/2+) 

0.047 

0.149 

0.962 

0.018 

0.650 

0.750 

-0.177 

0.977 

-0.115 

-0.279 

0.950 

-0.152 

Ac(l/2+) 

0.999 

0.001 

-0.020 

0.999 

<0.001 

-0.012 

0.917 

0.137 

0.374 

0.877 

0.289 

0.382 

Ac(1/2+) 

0.017 

0.100 

0.993 

0.010 

0.361 

0.931 

-0.138 

0.989 

-0.059 

-0.257 

0.957 

-0.132 

A6(1/2+) 

0.999 

<0.000 

-0.003 

0.999 

<0.001 

-0.004 

0.915 

0.141 

0.378 

0.876 

0.286 

0.390 


Significant mixing occurs only in the 1 / 2 + sector, for all flavors, particularly in the Sturmian models. Table lTvl 
shows the wave function coefficients for the two lowest 1 / 2 + states, in each flavor sector, for all four models 
(in the case of the nucleon, we show only the results from the HO models). The mixing shown in this table 
complicates the extraction of the form factors. However, in all results that we show for the form factors and 
the decay rates, this mixing is properly accounted for. Note that in each of these wave functions, there is also 
some contribution from the term in 774 . However, this component of the wave function has negligible overlap 
with the wave function of the parent baryon, and so is neglected here. 


B. Form Factors and Decay Rates 

In our calculation of the form factors, we have assumed that we can use non-relativistic approximations for 
the operators. This means that we have ignored terms in the various quark model operators that appear at 
order l/m^, l/mg, and above. Such terms have also been ignored in writing the hadronic matrix elements. 
However, in extracting the form factors, we have kept, and shown, terms that are of order l/(mqmQ). To 
examine the validity of this treatment, we write each form factor as 

F, = + + —FP+ 

rriq TUQ rngiTiQ 

= Ff ^ (61) 

and show the values for etc., in Table M In this table, we show only the results for the HONR and 

STNR models. 


For the elastic decays, the form factors Fi and Gi are dominant, while all other form factors are sub dominant. 
For 1/2“ final states, F 2 , Gi and G 2 are dominant, while for 3/2“, Fi and Gi are the dominant form factors. 
In each case, we see that the or term is significantly larger than the ‘higher order’ terms, as expected. 
The numbers in this table suggest that the convergence in l/iriq is rapid, modulo the model dependence. 
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TABLE V: Form factor components Ti and Qi as defined in Eq. EU, evaluated at the non-recoil point. The components 
are shown for the HONR (HO) and STNR (St.) models. The columns labeled ‘Ac’ are for the Ac —> A^*^ form factors, 
while those labeled ‘A^’ are for the At —> A^*^ form factors. 


form 

factor 

H.O. 

II 

< 

1/2+ 

H.O. St. 

H.O. 

> 

II 

1/2- 

Afo 

H.O. St. 

J+ = 
Ac 

H.O. St. 

3/2- 

A6 

H.O. St. 


0.98 

0.97 

0.99 

0.99 

0 

0 

0 

0 

-1.08 

-1.48 

-1.16 

-1.38 

T'f' 

0.54 

0.78 

0.20 

0.28 

0.36 

0.32 

0.16 

0.12 

-0.46 

-0.76 

-0.23 

-0.25 

^(Q) 

0.23 

0.15 

0.08 

0.05 

-0.04 

-0.04 

-0.04 

-0.01 

-0.10 

-0.37 

-0.03 

-0.12 

^(gQ) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-r(O) 

'^2 

0 

0 

0 

0 

-1.24 

-1.71 

-1.34 

-1.60 

0 

0 

0 

0 

'^2 

-0.54 

-0.72 

0.20 

-0.26 

0.36 

0.32 

0.16 

0.12 

0.46 

0.76 

0.23 

0.25 

'' 2 

0 

0 

0 

0 

-0.34 

-0.43 

-0.11 

-0.14 

<0.01 

<0.01 

<0.01 

<0.01 

'^2 

0.05 

-0.03 

0.01 

<0.01 

0.06 

0.05 

0.01 

0.02 

0.08 

0.07 

0.02 

0.01 

/p(O) 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

'r(g) 

3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

'^(Q) 

3 

-0.21 

-0.11 

-0.07 

-0.04 

0.34 

0.43 

0.08 

0.14 

0.37 

0.43 

0.13 

0.15 

/r*(gQ) 

'' 3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-rCQ) 

- 

- 

- 

- 

- 

- 

- 

- 

-0.14 

-0.13 

-0.06 

-0.05 


0.98 

0.97 

0.99 

0.99 

1.24 

1.71 

1.34 

1.60 

-1.08 

-1.48 

-1.16 

-1.38 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

g(.Q) 

0 

0 

0 

0 

0.04 

0.02 

0.02 

0.01 

0.07 

0.06 

0.03 

0.02 

gW) 

0.02 

-0.01 

<0.01 

<0.01 

0.06 

0.02 

0.01 

<0.01 

0.05 

0.07 

0.01 

0.01 

w~ 

0 

0 

0 

0 

-1.24 

-1.71 

-1.34 

-1.60 

0 

0 

0 

0 


-0.54 

-0.72 

-0.20 

-0.26 

0.36 

0.32 

0.16 

0.12 

0.46 

0.76 

0.23 

0.25 


0 

0 

0 

0 

0.08 

0.06 

0.04 

0.03 

0 

0 

0 

0 

gOlQ) 

-0.14 

0.06 

-0.02 

-0.01 

0 

<0.01 

<0.01 

<0.01 

0.14 

0.20 

0.02 

0.02 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0.23 

0.11 

0.08 

0.04 

0.08 

0.07 

0.04 

0.03 

-0.37 

-0.43 

-0.13 

-0.15 

g(qQ) 

0.13 

0.09 

0.02 

0.01 

-0.06 

-0.02 

-0.01 

0.01 

-0.17 

-0.1 

-0.03 

-0.02 

g(Q} 

- 

- 

- 

- 

- 

- 

- 

- 

0.14 

0.18 

0.06 

0.06 

gW) 

- 

- 

- 

- 

- 

- 

- 

- 

0.12 

0.13 

0.02 

0.02 


1 . Ac^A(*) 

In Table IVTl we show the values of the form factors at the non-recoil point, for the decays Ac —> A, for both 
elastic and inelastic channels. In this table, the results from all four models are presented. The results we obtain 
for the elastic channel are consistent with the predictions of HQET as estimated by Scora . 


In their treatment of the process Ac ^ Ae+J^, the CLEO Collaboration have used the leading order predictions 
of HQET to analyze the decay rate in terms of two form factors, and ^ 2 - In terms of the form factors that 
we have been using, these HQET form factors are 

a = + F2/2, 6 = ^^2/2, 

6 = Gi - G 2 / 2 , 6 = G 2/2 (62) 

The two sets of equations above arise from inverting Eqs. (EOl) either in terms of the E) or the G^. In Table IVTIl 
we show the values we obtain for the ratio evaluated at the non-recoil point. We also show the value 
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TABLE VI: The form factors for Ac 


a(*) 


transitions, calculated at the non-recoil point, in the four models used here. 


Spin 

model 

Fi F 2 Fs Fi Gi G 2 Gs Gi 

1 / 2 + 

HONR 

1.75 -0.54 -0.23 - 0.98 -0.54 0.23 - 

1 / 2 + 

HOSR 

1.76 -0.55 -0.24 - 0.98 -0.55 0.24 - 

1 / 2 + 

STNR 

1.90 -0.72 -0.11 - 0.97 -0.72 0.11 - 

1 / 2 + 

STSR 

1.78 -0.66 -0.09 - 0.92 -0.66 0.09 - 

1/2“ 

HONR 

0.32 -1.22 0.34 - 1.20 -0.80 0.08 - 

1/2“ 

HOSR 

0.42 -1.02 0.30 - 1.14 -0.61 0.10 - 

1/2“ 

STNR 

0.28 -1.82 0.43 - 1.73 -1.42 0.07 - 

1/2“ 

STSR 

0.36 -1.30 0.31 - 1.38 -1.04 0.08 - 

3/2“ 

HONR 

-1.83 0.46 0.37 -0.14 -1.00 0.46 -0.37 0.14 

3/2“ 

HOSR 

-1.81 0.52 0.35 -0.18 -0.94 0.52 -0.35 0.16 

3/2“ 

STNR 

-2.61 0.76 0.43 -0.13 -1.42 0.76 -0.47 0.13 

3/2“ 

STSR 

-2.03 0.58 0.34 -0.13 -1.11 0.57 -0.38 0.13 


obtained by the CLEO Collaboration in their analysis. We note that CLEO present a single value for the ratio 
of form factors, while we have two sets of values, arising from the two equations above. These two expressions 
give values for this ratio that are different, but not disturbingly so. The vector ratio (involving the Fi) tends 
to be smaller than the axial-vector ratio (involving the Gi), and both are smaller than the ratio extracted by 
the CLEO collaboration. The differences among the numbers we obtain using the two methods can be traced 
back to the l/mg terms in Ei; if those terms are ignored, both methods give the same value for the ratio. 


TABLE VIL The ratio C 2 /C 1 for ^(1/2^)- The first row is obtained using the vector relation defined in the text, 

while the second row is obtained using the axial-vector relation. 


6/^1 

HONR HOSR STNR STSR CLEO 

Vector 

-0.18 

-0.18 

-0.23 

-0.23 

-0.31 

Axial Vector 

-0.21 

-0.22 

-0.27 

-0.26 

-0.31 


Figure n shows the dependence of the form factors for the elastic transition Ac ^ A(l/2+), calculated in 
the HONR and HOSR models on the left, and in the STSR and STNR models on the right. In each panel, the 
solid curves arise from the SR version of the model, while the dashed curves are from the NR version. If we 
compare the form factors shown in Figure ^ we see that those calculated using the Sturmian wave functions 
have larger slopes near the non-recoil point (maximum than those calculated using the harmonic oscillator 
wave functions. The form factors calculated in the different models all have similar values near the non-recoil 
point (as seen in Table EH). The larger slopes in the case of the Sturmian model form factors means that we 
can expect smaller integrated rates from the STSR and STNR models. 

The differential decay rates, dT/dq^, that we obtain in the four models are shown in Figure [21 For these 
rates, we use \Vcs\ = 0.974. In these figures, we show the differential rates for decays to the elastic channel, as 
well as for two orbital excitations, the states with = 1/2“ and 3/2“. We have also examined the differential 
decay rates to the 3/2+ and 5/2+ orbitally excited states, as well as to the 1/2+ radially excited state. With 
the exception of the latter, we find these rates to be significantly smaller than those shown in this figure. 

As expected from the plots for the form factors, the differential decay rates that arise from the Sturmian 
wave functions for the ground state show a larger variation over the allowed q^ range. We also point out that 
the most noticeable difference between the NR and SR versions of a particular model is seen in the differential 
rate for the elastic decay. 

The integrated decay rate for the different final states in the different models are shown in Table EnH As 
anticipated above, the total semileptonic decay rates that we obtain in the harmonic oscillator models are 














Form Factors 


21 




[_^_,_I_^^^^_,_I_^^ 

0.0 0.3 0.6 0.9 1.2 

qVeV^) 


FIG. 1: Form factors for Ac ^ A(l/2^) obtained using harmonic oscillator wave functions (left panel, HOSR and HONR 
models) and Sturmian wave functions (right panel, STSR and STNR models). In each panel, the solid curves arise from 
the semirelativistic version of the model, while the dashed curves arise from the nonrelativistic version. Note that F 2 is 
indistinguishable from G 2 in all cases. 


TABLE VIII: Integrated decay rates for Ac —> A*-*^ in units of for different A states in the four models we 

consider. The last row shows the ‘elastic fraction’ obtained in our model, where the decays shown in the table are 
assumed to saturate the semileptonic decays. 


Spin 

r(HONR) r(HOSR) 

F(STNR) r(STSR) 

Expt. [311 

1/2+ 

2.10 

2.36 

0.79 

1.11 

1.05± 0.35 

1/2“ 

0.19 

0.29 

0.12 

0.15 

- 

3/2“ 

0.05 

0.06 

0.06 

0.05 

- 

1/2+ 

0.02 

0.02 

<0.01 

<0.01 

- 

total 

2.36 

2.73 

0.97 

1.31 

- 

TA/rtotal 

0.89 

0.86 

0.81 

0.85 

1.0 (assumed) 


significantly larger than those obtained in the Sturmian models. This effect is largest in the elastic decays, 
where the HO models predict decay rates that are more than twice as large as the ST models. We note that 
the elastic rates predicted by the ST models are much closer to the experimentally reported rate |3^ than those 
predicted by the HO models. 


From Table IVIIII it is clear that, while the elastic channel dominates the decay rate of the Ac, it does not 
saturate the decay. In each model, we find that the decay rate to the 1/2“ state is roughly one tenth of the 
elastic decay rate, while the decay rate to the 3/2“ state is about five percent of the elastic. Decays to these 
two excited states account for about 15% of the total decays of the Ac, assuming that decays to other excited 
states are negligible. It is also interesting to note that the ratio TA/Ttotai is almost independent of the model 
that we use, even though the absolute rates are very different in the different models. 

The assumption that the channels we explore saturate the resonant decays of the Ac is certainly consistent 
with the results we have obtained with the other states that we consider. First we point out that phase space 
limits how many excited A states can be considered, and the higher the excitation, the more limited the phase 
space available for producing such a state. For some final states for which there might be sufficient phase 
space to allow the decay, the spin-space structure of the state allows little overlap with the initial baryon, and 
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FIG. 2: The differential decay rates for different Ac ^ transitions, in the different models that we use. The curves 
on the left arise from the two versions of the harmonic oscillator model, while those on the right are from the Sturmian 
models. The curves are for exclusive final states with = 1/2”'’, 1/2“ and 3/2“. Also shown are the differential decay 
rates obtained by adding the exclusive modes described (labeled as ‘total’). In each panel, the solid curves arise from 
the semirelativistic versions of the models, while the dashed curves arise from the nonrelativistic versions. 


configuration mixing that could involve components with larger overlap with the initial baryon is very small. 
In addition, angular momentum factors (in orbitally excited states) lead to suppression of the decay rate. 

We can compare our predictions for decays to the excited A states with the assumption made by the CLEO 
Collaboration [^, that the elastic channel saturates the semileptonic decays of the A^. In our models, we find 
that between 11% and 19% of the A^ semileptonic decays are to excited states. In addition, our branching 
fraction (of 81% to 89%) to the ground state A must represent an upper limit, as we have not included any 
non-resonant production of multi-particle final states. It appears difficult to understand the lack of evidence 
for any decays to excited states in Ref. [^. This article reports no signal for decays of the kind Ac —> AXe~^v, 
and this is taken as evidence of saturation. However, the excited A states that we consider do not decay to Att, 
the most obvious decay mode to search for, as this decay is isospin violating. They will predominantly decay to 
Ett final states. In fact, the 1/2“ state, the A(1405), has a 100% branching fraction to Ett, while the A(1520), 
the 3/2“ state, has roughly equal dominant branching ratios to Ett and NK, with only about ten percent going 
into Atttt. Thus, our suggestion is that CLEO should investigate final states like ’ETriiy and NKiv^ and not 
states like Amrii'. 

The results discussed above are obtained using the assumption that the lightest of the = 1/2“ A states, 
identified with the S'oi state A(1405) found in analyses of scattering data, is a three-quark state. There are a 
number of other descriptions of this state in the literature, such as a dynamically generated bound state |45|. 
and a multi-quark state . If the CLEO Collaboration (or other groups) search for decays of the Ac to excited 
A states, especially the A(1405), and find no such decays, this would be a strong hint that this state is not a 
simple three quark state, as we have assumed. 

Our estimate of the fraction of Ac decays to excited states has important consequences for the absolute 
normalization of the branching fractions to the more than sixty observed final states in A+ decay. Most of these 
branching fractions are measured relative to the decay mode A/" —> pK~'K~^, and the absolute branching fraction 
of this mode cannot be extracted from data without introducing model dependence. One of the two important 
techniques for this extraction is based on measurements cross section for A+X production in e“''e“ 

annihilation, with the subsequent semileptonic decay A/" ^ At^Vi. The extraction relies on the assumption 
that the fraction / of decays A+ —> Xsi'^vi that have Xg as the ground state A is unity (the elastic channel 
saturates the semileptonic decays), with a significant uncertainty. Our calculated value / = 0.85, with an error 
of 0.04 estimated by evaluating / in four different models, changes the central value of this parameter and may 
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TABLE IX: Form factors of At, ^ Ac*\ calculated at the non-recoil point, in the four models we use. Also shown are the 
lattice estimates for the elastic form factors, taken from Ea. The lattice numbers are in fact multiples of their estimate 
of C(w), for which they explore a number of scenarios. 


r 

model 

Fi 

F2 

Fs 

F4 

Gi G2 

Gs 

G4 

1/2+ 

HONR 

1.27 

-0.20 

-0.08 

- 

0.99 -0.20 

0.08 

- 

1/2+ 

HOSR 

1.24 

-0.18 

-0.08 

- 

0.97 -0.18 

0.08 

- 

1/2+ 

STNR 

1.28 

-0.26 

-0.04 

- 

0.98 -0.26 

0.04 

- 

1/2+ 

STSR 

1.20 

-0.22 

-0.03 

- 

0.92 -0.22 

0.03 

- 

1/2+ 

Lattice 

1.28T0.06 

-0.19±0.04 

-0.06l°;“ 

- 

0.99 -0.24:+°°l 

0.09T0.02 

- 

1/2- 

HONR 

0.12 

-1.20 

0.11 

- 

1.21 -1.05 

0.03 

- 

1/2- 

HOSR 

0.15 

-0.95 

0.09 

- 

1.01 -0.82 

0.04 

- 

1/2- 

STNR 

0.10 

-1.63 

0.14 

- 

1.61 -1.50 

0.03 

- 

1/2- 

STSR 

0.11 

-1.21 

0.10 

- 

1.24 -1.12 

0.03 

- 

3/2- 

HONR 

-1.33 

0.17 

0.13 

-0.06 

-1.03 0.17 

-0.13 

0.06 

3/2- 

HOSR 

-1.13 

0.15 

0.12 

-0.05 

-0.87 0.15 

-0.12 

0.05 

3/2- 

STNR 

-1.75 

0.25 

0.15 

-0.05 

-1.36 0.25 

-0.22 

0.05 

3/2- 

STSR 

-1.31 

0.16 

0.11 

-0.05 

-1.04 0.16 

-0.18 

0.05 


allow a reduction in the assumed error from model dependence in the extracted absolute branching fractions. 


g. At ^ Ai*> 

In Table IDCI we show the values of the form factors at the non-recoil point, for the decays Af, —s- Ai*\ where 
this notation means that the Ac may be in an excited state. The results from all four models are shown, along 
with the results from a lattice st udy Ea The lattice results are actually given as multiples of ^{w), evaluated 
at the non-recoil point, and Ref. reports a number of different values for ^{w). In the ‘physical’ limit, values 
= i.03;°i« and = 0.87 ± 0.22 are quoted, where the two extractions are from the axial and 

vector currents, respectively. The results we obtain for the elastic decays are consistent with the predictions of 
HQET as estimated by Scora [^, as well as with these lattice simulations. 


Figure 131 shows the dependence of the form factors for the elastic decay of the Af,, calculated in the HONR 
and HOSR models on the left, and in the STSR and STNR models on the right. In each panel, the solid curves 
arise from the SR version of the model, while the dashed curves are from the NR version. As we noted in the 
case of the Ac ^ A(l/2+), the form factors obtained in the Sturmian basis have significantly larger slopes than 
the corresponding form factors calculated in the harmonic oscillator basis, at the non-recoil point. 

In terms of the Isgur-Wise function ^{w) for the elastic decay of the Af,, the form factor Fi is 


Fi 


1 + A 




(63) 


where A = mA^ — mb = m\^ — me at leading order in the heavy quark expansion. From the forms given in 
Appendix [nj and with the identification A « 2ma, we can extract 

in the harmonic oscillator basis, or 


i{w) 


1 


1 + 


p2 


(65) 
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FIG. 3: Form factors for A;, —> Ac(l/2^) obtained using harmonic oscillator wave functions (left panel, HOSR and HONR 
models) and Sturmian wave functions (right panel, STSR and STSR models). In each panel, the solid curves arise from 
the semirelativistic version of the model, while the dashed curves arise from the nonrelativistic version. Note that F 2 is 
indistinguishable from G 2 in all cases. 


in the Sturmian basis (assuming single-component wave functions), and we have assumed that a\ = ax' = a, 
P\ = /3a' = /3 in the heavy quark limit. Writing 


- 1 ) « - 1 ), 

the above expressions become 

Ciw) « exp - 1)^ 


in the harmonic oscillator basis, or 


C(w) 


1 

1 I 3m2 / ^ N 

1 -F -0r{w — 1) 


in the Sturmian basis. 

The Isgur-Wise function may be expanded as 

^{w) = 1 - p^{w - 1) -F y (w - 1)^ -F ..., 


( 66 ) 


(67) 


( 68 ) 


(69) 


where the slope of the form factor at the non-recoil point has been denoted p^, and the curvature is denoted 
cr^. Rigorous bounds have been placed on the values of both the slope and curvature parameters for meson 
decays, and some models have difficulty in satisfying those bounds. In particular, in the model of ISGW 0, a 
factor K was introduced by hand (see the discussion between Eqs. (B2) and (B3) of Ref. Q) to modify the 
dependence of the form factors. In our model, the equivalent procedure would be to change Ih in Eq. m 
from 


' 3/2 3/2' 
ad ad, 


exp 


^AA' 


aly j 


Ih = 
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as calculated to 


Ih 


' 3/2 3/2' 

a 1 a 1 , 


exp 


'^AA' 


( 3m^ \ 

V j ■ 


The argument used by ISGW was that this factor of k would take into account ‘relativistic effects’. The effect 
of this change is shown in Figure^ where the form factors for A;, —> are plotted as functions oi w = v ■ v', 
for the two harmonic oscillator models (upper graphs). For comparison, the lower graph shows form factors 
obtained in the Sturmian basis, also as functions of w. The graph on the upper left shows our calculated form 
factors, while that on the upper right shows form factors including the factor of k. 





w 


FIG. 4: The elastic form factors for the decay of the A(, as functions of w. The upper two graphs arise from the harmonic 
oscillator model, while the lower graph is from the Sturmian version of our model. Among the upper graphs, the panel 
on the left shows the form factors obtained in this work, while those on the right incorporate a ‘relativistic’ factor in the 
exponential (see text). In each panel, the solid curves arise from the semirelativistic version of the model. 
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Table H lists the slope of the Isgur-Wise function that we have extracted, at the non-recoil point, in both the 
harmonic oscillator and Sturmian models, as well as in the ‘relativistically modified’ harmonic oscillator model 
(using the factor a). The slopes of the form factor near the non-recoil point are larger in the Sturmian models 
than in the harmonic oscillator models. This is easily understood by noting that the value of is 


in the HO models, and 




(70) 


(71) 


in the ST models. The extra factor of two in the latter case arises because the form factors in the ST models 
have a dipole dependence on w. A corresponding monopole form would give the same slope as the HO models. 
Since the values of nia- are similar in the two sets of models, and the values of a are not very different from the 
values of f3, the ST models will give slopes that are roughly twice as large as the HO models. In the same way, 
it is easily shown that the ST models lead to curvatures that are about six times as large as those obtained in 
the HO models. 


TABLE X: Slope of the Isgur-Wise function, evaluated at the non-recoil point, for the elastic decay of the A^. 


model 

HONR HOSR HONRk HOSRk STNR STSR 

d^(w)/dw 

-1.38 -1.33 -2.82 -2.71 -5.71 -3.27 


The K-modified harmonic oscillator model leads to slopes that are similar to those obtained in the Sturmian 
models, since the value chosen for k was 0.7 (so that 1 /k^ « 2). Relativistic effects do not need to be invoked 
to obtain the large slopes obtained in the Sturmian models. The differences in the slopes are simply artifacts 
of the expansion bases used for the wave functions. 

In the follow-up article to Ref. 0, Scora and Isgur rewrite the quark model form factors, explicitly 
replacing the exponential factor that arises with the harmonic oscillator wave functions. The change they make 
is 


exp 


-IrlflimB - muY 

D 



1 

{1 + ’ 


where is the value obtained from the harmonic oscillator wave functions, and 


(72) 


r 


2 


3 

AmQmq 


+ ^wf + ’’qcdj 


(73) 


where the last term arises from matching of currents in HQET with full QCD. In Eq. the integer N = 
2-\-n + n' ^ where n and n' are the harmonic oscillator principal quantum numbers for the initial and hnal wave 
functions. The hnal forms that they used are therefore very similar to the forms that we have obtained in the 
Sturmian models. 

The values we have obtained for the slope of the Isgur-Wise function in our Sturmian models are signihcantly 
larger than the value obtained recently by Huang et al. using a HQET approach based on QCD sum rules: 
their value for is less than 1.5, similar to the values we obtain in the HO models. In a recent analysis of 
the Ab form factor measured in hadronic Z decays, the DELPHI Collaboration found p^ = 2.03 ± 0.46, 
where the error shown is statistical. They also reported two sets of systematic errors, each comparable to the 
statistical error. This result means that for the Sturmian models, we will obtain integrated decay rates that 
are significantly smaller than the DELPHI rate. In the lattice study by Bowler et al. the reported slopes 
is 1.1 ± 1.0. A more recent lattice study with O (a^, asa^) improved lattices 0 does not quote values for the 
slope. However, a conservative estimate from the graphs they present gives values for that appear to be 
consistent with the large values we obtain in the ST models. 







dr/dq (s GeV ) dr/dq (s GeV ) 
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Also of some interest is the curvature of the Isgur-Wise function, denoted cr^. In the HO models with no 
modifications, the prediction is that cr^o = (Pho)^i while the ST models give cr|q^ = 3(p3qn)^/2. Bounds on the 
curvature of the Isgur-Wise function for meson decays have been derived by Le Yaouanc, Oliver and Raynal [^. 
To the best of our knowledge, no such bounds have been derived for baryon decays. However, the values of 
the curvature we obtain using both the HO and ST models easily satisfy the known bounds for meson decays. 
Note that the large slope and large curvature we obtain suggest that the common procedure of parameterizing 
the Isgur-Wise function only in terms of its slope parameter, can potentially lead to significant errors in the 
extraction of CKM matrix elements. 




q^GeV^ q^GeV^ 




q^GeV^) q^GeV^) 


FIG. 5: The differential decay rates for different Ab —> Ac*^ transitions, in the various models that we use. The curves 
on the left arise from the two versions of the harmonic oscillator model, while those on the right are from the Sturmian 
models. The upper panels are for Ab —> Ai*'^iue, where i is e~ or /i“. The lower panels are for Af, —> A[^\vt- The curves 
are for final states with = 1/2’'", 1/2“ and 3/2“. 

The differential decay rates dT/dq^ that we obtain in the four models are shown in Figure 0 (assuming 
\Vcb\ = .041). In these plots, we show the differential rates for the elastic channel, for the radially excited 1/2)'’ 
state, as well as for decays to two orbital excitations, the states with = 1/2“ and 3/2“. We have also 
examined the decay rates to the 3/2"'", 5/2"'’ states, and found them to be smaller than those shown in this 
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figure, contributing of the order of one or two percent to the total rate. 

The integrated decay rate for the different final states in the different models are shown in Table EH As 
anticipated above, the total semileptonic decay rates that we obtain in the harmonic oscillator models are 
significantly larger than those obtained in the Sturmian models. This effect is largest in the elastic decays, 
where HO models predict decay rates that are more than twice as large as the ST models. Note that, in all 
models, the decay rate to the 3/2“ state is roughly twice the decay rate to the 1/2“ state. In the heavy quark 
limit, this ratio of decay rates is expected to be two, and results from arguments that are similar to spin-counting 
arguments. 


TABLE XL Rates for At —> decays in units of The first five rows are for decays with a muon or electron 

in the final state, while the last four rows are for decays with a r in the final state. The rows labeled ‘total’ are obtained 
by adding the exclusive decay rates shown in the table, while the row with the branching fractions assumes that the 
exclusive channels shown saturate the semileptonic decays of the At,. The elastic fraction reported by the DELPHI 
collaboration (fifth row, sixth column) is actually r(A|,-»A ?i^)+r(At^A ■wTreni) ’ errors on both DELPHI results are 
statistical and systematic, respectively. 



r(HONR) r(HOSR) r(STNR) r(STSR) 

Edelphi 

1/2+ 

4.60 

5.39 

1.47 

2.00 

4 r»'^+0.90+1.30 
-0.65-0.98 

1/2“ 

0.45 

0.52 

0.26 

0.27 

- 

3/2“ 

0.95 

0.91 

0.63 

0.61 

- 

Total {Ai*'^£~i'i) 

5.95 

6.82 

2.36 

2.88 

- 

r Ac /r total 

0.76 

0.79 

0.62 

0.69 

r\ yj^^+0.10+0.07 
-0.08-0.06 

1/2+ 

1.90 

2.09 

0.82 

1.00 

- 

1/2“ 

0.10 

0.11 

0.08 

0.07 

- 

3/2“ 

0.15 

0.13 

0.14 

0.12 

- 

Total (A/*%“!?t) 

2.15 

2.33 

1.04 

1.19 

- 


Table Ixll also shows that a significant fraction of the semileptonic decay of the At is inelastic. This is analogous 
to what has been seen in B semileptonic decays, where the elastic channels account for no more than about 
80% of the total semileptonic decay rate. For the At, our predicted ratios are similar, ranging from 62% to 77% 
of the total semileptonic decay rate. We have estimated the total semileptonic decay rate by assuming that the 
three exclusive modes shown in Table ED saturate the semileptonic decays (rates to other states that we have 
examined are significantly smaller than those shown in the table). Using these numbers, we obtain predictions 
for the total semileptonic decay rate of the At, also shown in Table ED 

For comparison, the PDG gives a rate of 7.486 ± 2.105 x 10^°s“^ for the inclusive semileptonic decay 
At —*■ Ac^v-\- anything. This is signihcantly larger than any of the total semileptonic widths we obtain, but the 
authors of the PDG emphasize that this value results from assumptions about the fragmentation of b quarks into 
baryons, and ‘cannot be thought of as measurements’ El The DELPHI value for the elastic semileptonic decay 
rate is also shown in Table ED As anticipated, the rates we obtain in the Sturmian models are signihcantly 
smaller than the DELPHI rate, while those obtained in the harmonic oscillator models are consistent with the 
DELPHI measurement. 

The above examination of the decays of the Ac found that the Sturmian models provided rates that were 
consistent with the CLEG measurements, while the harmonic oscillator models gave rates that were twice as 
large. This suggested that the Sturmian models might be more reliable. For the A;, decays, we see that the 
harmonic oscillator models provide rates that are more consistent with the single measurement available to date. 
For the Sturmian models, the predicted rates are about 2a away from the reported value, if the systematic and 
statistical errors are treated in quadrature. 

The DELPHI Collaboration also reported on the elastic fraction of the semileptonic decays of the Af,. For the 
ratio Y{Ai^Ace^e)+'^{'th^KcT^Trtve) ’ ^ value of 0.47lQ;JglQ;Qg, with no evidence for resonant decays. This 

ratio is smaller than we predict, in all models. However, our predictions must be thought of as upper limits 
for the elastic fraction, as we do not include any non-resonant semileptonic decays. We note that our predicted 
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ratios are already somewhat smaller than those reported in the decays of B mesons, while the DELPHI ratio 
is smaller still, suggesting that there are signihcant differences between the semileptonic decays of the heavy 
baryons and those of the heavy mesons. If the DELPHI results for both for the elastic rate and the elastic 
fraction are not modified by future experiments, this aspect of the physics of heavy hadrons will require further 
scrutiny. 


3. Aq —> Decay 

The decays of the Aq to final states consisting solely of light quarks are interesting as they provide an alternate 
means of extracting CKM matrix elements like Vub- The expectation from HQET is, modulo 1 /toq effects, that 
the form factors that describe the Ac —> n semileptonic decays will be the same as those describing the At, ^ p 
semileptonic decays. To explore this, we now examine the form factors for these two decays. 

In Pigure|n|we show the form factors and ^2 for the transitions Ac —*■ n and At, ^ p, obtained in the 

two harmonic oscillator models. The two forms are found using the two sets of equations in Eq. 16211 . The 

value of ^2 is independent of which of the two sets of equations we use, up to the order to which we calculate 
the form factors. In both the nonrelativistic and semirelativistic versions of the model, the two curves for 
(top right plot in Fig. are very similar, indicating that the HQET prediction, that this form factor should 
be the same for both transitions, indeed holds up to small corrections. For the semirelativistic version, the two 
curves are closer than in the nonrelativistic case. The differences seen in the curves for ^ 2 , which are consistent 
with those in the curves for arise mainly from the differences in the size parameters (ap and av) between 

the Ac and At, states in the models (see Table I hI. The curves for (top left plot in Fig.EJ show the biggest 
differences in going from A;, —> p to Ac ^ n, in both models. Here, the differences get some contribution from 
the 1 /toq term that is present in Fi. 

In Figure [71 we show the differential decay rates for Af, decaying semileptonically into the four lowest-lying 
nucleon states, while Table Eul shows the integrated rates into six exclusive states. Also shown in this figure 
and table are the rates that we obtain when the final lepton is a r. The ground state nucleon is the largest of 
the CKM suppressed decays of the At,, but it accounts for less than 50% of these decays, in both of the harmonic 
oscillator models. A large fraction (about 20%) goes into the first excited state, the Roper resonance, usually 
treated as a radial excitation of the ground-state nucleon, as it is in this model. As with the A(1405) in the 
decays of the Ac, this result hinges on the assumption that the Roper resonance is a three-quark state, and that 
it is the first radial excitation of the nucleon. A number of hypotheses for the internal structure of this state 
have been made, such as pentaquark partner [^, dynamically generated state |^, and hybrid state |5^. In 
each of these scenarios, the rate at which the At, decays semileptonically into this state is affected by its internal 
structure. For the three-quark, radially-excited scenario, the prediction is that decays to this state are about 
60% of the decays to the ground state nucleon, a rather large fraction. If ample A;,’s can be produced, their 
semileptonic decays may therefore provide information that can be used in understanding the structure of the 
Roper resonance. 

We have examined decays to other excited nucleons, and those shown in Table ixm are by far the dominant 
ones. We have also examined one additional 1/2+ nucleon state, two additional nucleon states with J+ = 3/2+, 
and one additional nucleon state with = 5/2+, none of which are shown in Table jXHl Of these, the rate to 
the additional 1/2+ state is less than 1% of the ‘total’ rate that we have estimated, while rates to the additional 
3/2+ and 5/2+ states are similarly small or even smaller. These small rates are a direct consequence of the 
structure of these states, as their overlaps with the decaying A^, in the spectator assumption, are very small. 
The only other excited nucleons that may occur with ‘significant’ rate in the semileptonic decays of the A;, are 
those with higher spins, such as 7/2+ and 5/2“. However, for such states, orbital angular momentum centrifugal 
factors will lead to some suppression of the decay rate. 


Figure 0 shows the differential decay rate for Ac ^ n, while the integrated decay rates for two exclusive 
modes Ac ^ obtained using \Vcd\ = 0.224, are shown in Table IXTTl It is clear from this table that decays 

of the Ac to excited states of the nucleon are strongly suppressed, due in part to the reduced phase space. 
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FIG. 6: Form factors (top left), (top right) and ^2 (bottom) for the transitions Ac —> n and At —> p. All curves 
are found using the harmonic oscillator models, with the solid curves corresponding to HOSR, and the dashed curves to 
HONR. The two plots for arise from the two ways of evaluating this form factor, shown in Eq. (HU. 


VII. CONCLUSIONS AND OUTLOOK 

A constituent quark model calculation of semileptonic decays of A;, and Ac baryons, which has several novel 
features, is described here. Analytic results for the form factors for the decays to = 1/2+ ground states and 
excited states with different quantum numbers are evaluated, and compared to HQET predictions. For A;, —> Ac 
transitions, the relations among the form factors, predicted by HQET, are satisfied by the form factors obtained 
in the model, independent of the basis used to describe the baryon wave functions. For the elastic form factors, 
as well as for the form factors for decays to the (l/2“, 3/2“) doublet, the HQET relationships among the form 
factors are found to hold up to the order we have examined, namely l/rrih and Ijmc- For states of higher spin, 
we have compared our model form factors to the HQET predictions at leading order, and the expected relations 
hold at that order. 






dr/dq (s' GeV ) dP/dq (s' GeV ) 
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FIG. 7: Differential decay rates for Ab in the HONR and HOSR models. The upper panels show the rates for 

Aj, —> e~Ve, while the lower panels show the rates for Ab t~V r^ both in units of |K,6|^. The panels on the 

left are from the HOSR model, while those on the right are from the HONR model. 


These form factors depend on the size parameters of the initial and final baryon wave functions, and so a fit 
to the spectrum of the states treated here is performed. Two model Hamiltonians are used, with either a non- 
relativistic or semi-relativistic kinetic energy term, and with Coulomb and spin-spin contact interactions. The 
wave functions are expanded in either a harmonic oscillator or Sturmian basis, up to second-order polynomials, 
and our numerical results for form factors and rates are calculated using the resulting mixed wave functions. 
Four sets of predictions are made for form factors and rates, with wave functions, size parameters and mixing 
coefficients arising from fits using both the non-relativistic and semi-relativistic Hamiltonians, and using the 
two different bases. These predictions can be used to assess the model dependence in the results we obtain. 

Interestingly, the form factors for decays to ground state daughter baryons evaluated using the Sturmian basis 
for the wave functions have slopes at the non-recoil point that are significantly larger than those evaluated using 
the harmonic oscillator basis. As a result, the corresponding integrated decay rates for Ac ^ A elastic decays. 
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TABLE XII: Decay rates of Ab —> in units of lO'^^s ^ x Also shown are the rates for Ac ^ N^*''°£'^ve, 

in units of obtained using \Vcd\ = 0.224. 



Ab ^ 

N<'*'>+rue 

Ab ^ 

N<'*'> + T-p-r 


r(HONR) r(HOSR) 

r(HONR) r(HOSR) 

1/2+ 

4.55 

7.55 

4.01 

6.55 

1/2+ 

2.92 

4.44 

2.20 

3.05 

1/2- 

1.42 

3.85 

1.10 

2.73 

3/2" 

1.54 

1.82 

1.03 

1.07 

3/2+ 

1.03 

2.16 

0.28 

0.58 

5/2+ 

0.79 

1.49 

0.38 

0.55 

Total 

12.25 

21.31 

9.00 

15.53 


Ac ^ 


- 

- 

1/2+ 

1.02 

1.35 

- 

- 

1/2- 

0.02 

0.04 

- 

- 



FIG. 8: Differential decay rate for Ac —> n in the HONR and HOSR models. 


calculated using the Sturmian wave functions, are smaller than those obtained using the harmonic oscillator 
basis wave functions. The Sturmian rates are both consistent within errors with the experimentally reported 
rate of 1.05 ± 0.35 x 10^^ s“^, while those calculated using the harmonic oscillator basis are significantly larger. 
As pointed out by Keister and Polyzou although calculations using the Sturmian basis are not as simple 
as those using the harmonic oscillator basis, the resulting form factors have shapes which are expansions in 
inverse powers of 1 + with k the decay three-momentum (in a non-relativistic decay calculation like 

ours), and A a constant which is calculated in terms of quark masses and wave function size parameters. This is 
closer to the form expected from experimental studies of hadron decay form factors, and so the use of Sturmian 
basis functions produces realistic results for decay calculations even with the inevitable truncations of the basis 
required for tractability. Larger scale numerical calculations using the Sturmian basis require fewer basis states 
than those using the harmonic oscillator basis to yield accurate energies and decay form factors for excited 
states. 

Although the use of a semi-relativistic Hamiltonian does not necessarily lead to a better fit to the spectrum, 
in calculations using both bases it results in an integrated decay rate for Ac ^ A elastic channel that is closer 
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to the central value that has been experimentally reported. However, the rate obtained in the nonrelativistic 
version of the Sturmian model is also consistent (within 1 a) with the experimentally reported value. 

Decay form factors and rates to all available excited state daughter baryons are evaluated using these four 
models. Significant branching fractions are found for Ac inelastic semileptonic decays in all four calculations, 
with the total to all excited states ranging from 11 to 19%. This has important consequences for the absolute 
normalization of the branching fractions to the many observed final states in A+ decay, most of which are 
measured relative to the decay mode A+ ^ pK~'K'^. The extraction of the absolute branching fraction of this 
mode, from measurements of the cross section for A+Ai production in e“'"e“ annihilation, requires knowledge 
of the fraction / of semi-leptonic decays A+ ^ to the elastic channel. This contradicts the available 

CLEO analyses, in which it is assumed that the elastic decay of the Ac saturates its semileptonic decays. A 
larger fraction, from 23 to 38%, of Af, ^ Ac semileptonic decays, are found to be inelastic. Elastic decays of the 
Ah involving tau leptons in the final state are suppressed by roughly a factor of two because of the reduction in 
the final-state phase space, and those to excited baryon states are suppressed more strongly. 

HQET predicts that the form factors and ^ 2 , defined earlier, should be the same for the decays A;, ^ N'^ 
and Ac ^ N~^, up to 1/mb or l/ruc corrections. Within our models, we find that the two form factors are very 
similar, but not identical, with the differences arising from differences in the size parameters for the Af, and Ac. 
In the case of the decay of the Af, to nucleons, we find that the ‘elastic’ fraction is quite small, of the order of 
35% when the leptons produced in the decay are light. A number of excited nucleons contribute to the total 
rate, with the radially excited Roper having the next largest branching fraction. This may be used as a test of 
the structure of this resonance, if ample Af,’s can be produced. 

The work presented in this manuscript can be extended in a number of directions. We plan to examine 
the semileptonic decays of heavy flq baryons, both to ground states and to a number of excited states, in a 
calculation similar to the one outlined here. Since the description of these states using both the quark model 
and HQET is more complex, it will be interesting to see if the correspondence between quark-model results 
and the predictions of HQET still holds. We can also apply our model to the description of the semileptonic 
decays of the light baryons, although these are already successfully described by Cabbibo theory. Essentially all 
experimentally accessible observables for these decays have been measured, and it will be interesting to see if 
our model, constructed with no special reference to chiral symmetry or current algebra, can describe the results 
of these measurements. 

We have not examined the predictions of our model for the many polarization observables which can, in 
principle, be measured in semileptonic decays. One example is the asymmetry parameter in the decays of 
the Ac, which has already been extracted by the CLEO collaboration. The predictions of our model for this and 
similar quantities are therefore of some interest. In addition, the rare decays of heavy baryons, such as Af, ^ A 
can easily be treated in the framework that we have developed. Such processes, along with their meson analogs, 
are used in searches for physics beyond the standard model. However, the interpretation of the measured rates 
depend strongly on estimates of the form factors involved (in much the same way that extraction of CKM matrix 
elements depend on the form factors that describe semileptonic decays). Finally, if factorization, in some form, 
is valid, the semileptonic form factors calculated in the manuscript may also be useful in the description of 
nonleptonic weak decays. 

It may also be possible to systematically improve the quark model used in the present calculation. An obvious 
first step is the implementation of full symmetrization of the spatial wave functions in the Sturmian basis, which 
would allow calculation of results for decays to final state nucleons in this basis. 
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APPENDIX A: WAVE FUNCTIONS 

As mentioned in the text, our wave functions are expanded in two different bases. For the states of different 
spins and parities considered here, the expansions are given in this Appendix. For Aq states with = 1/2“'", 
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the expansion reads 


'^Kq,i/ 2 +m — '?^Aq ’7f'^V'oooooo(Pp,P a) + %®'0ooiooo(Pp,P a) + ?7^®V’ooooio(Pp,P a) Xi/ 2 i^) 
+ ?74®V'00010i(Pp,Pa)Xi/ 2W + [V'lMi,010l(Pp,PA)xf/2(^ - ^i)] 
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where [V'LMpnp^prix^>(PpjP a)xs(-^ ~ jm ^ shorthand notation that denotes the Clebsch-Gordan sum 
- ML)lpLMLnpipnxex{Pp,P\)Xs{M - Ml)- 
For Aq states with = 1/2“ and 3/2“, the expansion is 
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where J can take the values 1/2 or 3/2. 

For Aq states with = 3/2+, the expansion is 
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J 3/2,M 
- 3/2,M 

(A3) 


For = 5/2+, the expansion is 

^Aq,5/2+M = </Aq (?7f‘^V'lMz,010l(Pp,PA)xf/2W + [V'2Mz,0101 (Pp, PA)xf/2 “ ^a) 


5/2,M 


Aq 

% 


Aq 

% 


■02Mi,OlOl(Pp, Pa)Xi/2(^ - ^a) 
V’2Mi,0002(Pp, Pa)x5'/2(^ - ^a) 


J 5/2,M 


Aq 

•^4 


V'2Mz,0200(Pp,Pa)x)’/2(^ - ^a) 


J 5/2,M 


J 5/2,M 

(A4) 


No other states are expected to have significant overlap with the decaying ground-state Aq in the spectator 
approximation that we use. 

The wave function components for nucleons are different from those shown above, due to the different (12) 
symmetry in the wave functions, and are shown below. For = 1/2+, nucleon wave functions are expanded 
as 


'^N,l/2+M = (^N [’7rV'000000(Pp,PA) + ??^V’ooiooo(Pp,Pa) + ?7:fV'00001o(Pp,PA)] Xl/2(-^) 

+ lA'4’mmm{Pp,Px)Xi/2i^)+lb V'imz,oioi(Pp,Pa)x)’/2(^ - ^a) 

[V'2Mz,0200(Pp, Pa)X3/2(-^ “ -^a)] ^ [V'2Mj:,0002(Pp, Pa)X3/2(^ “ ^a) 


1q 


1/2,M 
.S 


(A5) 


J 1/2,M 
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For = 1/2 and 3/2 , the expansion is 




N,J-M 


(t^N ( 


+ri2 


V3 


Vi 


J JM 


V'iMiOioo(Pp,PA)Xi/2(M - Ml) 
■f/lMpODOllPp, PA)xf/2(^ - Ml) 
ipiM,^Q00iiPp,P\)Xi/2iM - Ml) 


JM 


where J can take the values 1/2 or 3/2. 


(A6) 


APPENDIX B: INTEGRALS IN THE STURMIAN BASIS 

Wave functions expanded in the Sturmian basis have been used by other authors in exploring aspects of heavy 
meson phenomenology |55|| . However, to the best of our knowledge, there have been no prior applications to 
baryon phenomenology. We therefore believe that it is useful to outline some of the steps needed in using this 
basis for calculations of the kind that we present. 


1. Integrals for Hamiltonian Matrix Elements 

We begin by reminding the reader that, in coordinate space, say, the spatial wave function components are 
written as 


l/iMriptpnxG (p, A) = '^{LM\£pm,exM - m)lpn,e,mip)'4’nxixM-mW, 

m 

with p and A as defined in the main text. 

In the Sturmian basis, evaluation of the matrix elements of the non-relativistic kinetic energy operator, as 
well as those of the parts of the potential that depend only on ri 2 = |ri — r 2 |, are relatively straightforward, 
in the latter case because p = However, the evaluation of terms that depend on ria or r 2 z is not 

as straightforward. To illustrate the way in which such calculations are carried out, we consider the linear 
potential, and examine the term 


P/“ = 5|ri-r3| = 6ri3. 


(Bl) 


We begin by writing 


D3 = 


V2 


p + V3A = IJ + A'l = + 2V3p • A + 3A2)1/2. 

v2 


In the above, p' = plx/2 and = a/ 3/2A. The latter form is expanded in spherical harmonics, yielding 


D3 


^ 0 J- 


i/2 


^ (2Z + 1) A"+i V(2^ + 3) (2Z-1) 

for p' < A/ and a similar expression with p' X' otherwise. In this expansion, 

(Yiip) ■ y,(A)) ^ ^(-i)-yr(p)A/-'"(A) 


(Yi{p)-YiCX)) 


(B2) 


(B3) 


(B4) 


Calculation of (ri 3 ) then requires the evaluation of the matrix element {L'n'pl'pn'^l'^\yi{p) ■ yi{X)\Lnplpn\lx), 
which symbolically denotes integrations over the angles defining p and A. This is done with the use of 6-J 
symbols, leaving integrals over the magnitudes of p and A which can be done either numerically or analytically. 
For the potentials we use, all terms can be handled analytically. Terms in the potential that depend on r 23 are 
handled in a similar manner. 
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2. Integrals for Current Matrix Elements 


In order to evaluate the form factors in the Sturmian basis, integrals of the form 

(P) yi-2 (P + ak) (p) 




_|_ q, 2^ (^P _|_ 


(B5) 


must be calculated. In the above, p represents an internal momentum conjugate to one of the Jacobi coordinates 
(for these integrals, px), while k is the momentum of the daughter baryon in the frame in which the parent is at 
rest. The constant a = with mjy^ being the mass of the daughter baryon in the decay. The 

quantities yi (p) are the vector harmonics, with £ 1^2 being the orbital angular momentum in the initial or final 
state, respectively, while y^^ (p) arises from the Pauli reduction of the vector or axial current. For simplicity 
we choose l\= = 0, but this will still be sufficient to illustrate the method. 

With the use of Feynman parametrization, this integral is first rewritten as 


7 - 0,^,0 


1 r(ni + 712 ) 

VItt r(ni)r(n2) 

1 r(ni + 712 ) 
VItt r(ni)r(n 2 ) 
1 F(ni + 712) 
VItt r(ni)r(n2) 


d.J 

f d?p 


f d?p 


f d^p 


(p) (p + gk) 

I a: + (1 — x) (p + ak)^ + | 

_ - x)^^-^yo (p) ye (p + gk) _ 

+ 2a(l — a:)p • k + — a;) + /^^(l — x) + 

^ni-l(l _ (-p^ (-p 


|[p + a(l — a;)k]^ + a'^k'^x{l — x) + + /^^(l — x)| 


ni+n2 


(B6) 


where the factor of I/VItt arises from one of the vector harmonics with £ = 0. 
Defining 


u = p + a(l — x)k 


(B7) 


and substituting into the integral gives 


jQX,o ^ 1 r(ni+n2) f 3 

V4^r(ni)F(n2)io i 


x”i ^(1 —x)”2 ^3^0 (u) 3^7 (u + axk) 


[■u^ + a?k"^x{l — x) + a^x + — x)]"^'*’"^ 

The angular integration can be performed after expanding the 34 (u + axk) to give 

x”i-i+^(l - x)”=^-i 


Ti 




dui 


+ a^fc^x(l — x) + a^x + /3^(1 — x)] 


ni+n2 ■ 


(B8) 


(B9) 


Using 


,2m 


du 


1 r(TO+l/2)r(n-m-1/2) 


[v? + Af 2 yI"-™-V 2 


r(n) 


(BIO) 


in the above equation gives 


n-i n N r(ni + 712) r( 3 / 2 )r(ni + 712- 3 / 2 ) 

7n,„ - “»(k) r(„,)r(„,) r(n. + „.) 


„ni — 1+^ 


dx- 


(1 -x) 


n2 —1 


lo 2 [a^k‘^x{l — x) + a‘^x + /?^(1 — 
x"i-i+^(l - x)"=-^ 


r(ni)r(n 2 ) Jo 2[a‘^k‘^x{l — x)a^x 

This integral can now be written as a sum of terms with 


(Bll) 


yjm _ 
*-'n — 


dx 


0 (co + CiX + C 2 x 2 )"'''^/^ 


(B12) 
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where 


cq= 0^, Cl = a?k'^ + — 0^, C 2 = —a^k^. (B13) 

Each of these terms can be then be integrated analytically to give the required matrix element. 

This procedure works as long as 2n > m. When 2n < m, the last integration leads to logarithms. Such terms 
are expanded around fc = 0 before the form factors are extracted. 


APPENDIX C: EXPRESSIONS FOR THE FORM FACTORS 

The analytic expressions that we obtain for the form factors are shown in the following subsections. The results 
shown are valid for single-component wave functions. We separate the results obtained using the harmonic 
oscillator basis from those obtained using the Sturmian basis. 


1. Harmonic Oscillator Basis 

a. 1/2+ 


where 


^AA' 


Fi = Ih 


1+ — 


E, = -I. 


H 


“aa' 

mcr Qfx 


Too 


mq 

a\a\, 


mqaly ialyrngiriQi’ 


F3 = —Ih 
G i = Ih 
G2 = —Ih 


TTla at 


TTIQ aly 


2 ^,2 


1 - 


at a 


A'^A' 


12a^;^,m,mQj 

rUrr 0 ?y 


-b 


A'-'A' 


mn a 


AA' 


G 3 = I, 


H 


TTZ/y GL \ 


12 mqmQaj^y^, 

mlalal, 1 


1-b 


12 mi 


rriQa^y rriqmQa^y 


Ih=(^ 
\ a 


^XX' 


3/2 


exp 


3mz 


2 m\ 


a 


' — ^<7 is the mass of the light quark. 


a 


XX' 


XX' 


b. 1 / 2 + 


Fi 

F 2 

F 3 


= Ih 




(al-al)- 


nia 


3'^AA' 


= -Ih 


6 mqa\y 


[7a1 — 5ay) 


^(7ai-3ai,) + ^ 

rriq mg 


(7ai,-3ai) 


Too - 


4mQj ’ 


ajjna /_ 2 q 2\ 
Ih ^ ^ 3 Go:y - 3aA), 


bmga 


AA' 
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Gi 

G2 

G3 


= Ih 
= -Ih 
= -Ih 


(al-al) 


2^,2 


aXa 


X^X' 


2q;aa' 


where 


72a\ymqmQ 


{7al-3al,) 


Qmqa\y L 
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QmQa\y 


Ih = 
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-) 


QmQj 
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axay 
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2 
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3mi 


’■w 


c. 1 / 2 “ 


where 
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TUq 

2m„ 
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TOqJ 

a\ 
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= Ih 
= Ih 
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ax QmQ 
2 ma , ax 


6mqmQa‘ly 
, «A 
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H 


ax 

SniQ 


1 - 


SmgJ ’ 


—^( 3 al- 2 al,) 
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Ih = 


OfAttA' 


a 


XX' 


5/2 


exp 


3 mi 


2mi 


"XX' 
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Fi 
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F3 

Fi 

Gi 


= Ih 


irria 

ax 




a 


XX' 


mn 


mg 


= -I. 


H 


3rnz 


a. 


7)axa\,ma 


= -Ih 
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3 m‘i 


ax 


ax 


niQ aly 2mQ\ 
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ax 
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where 


where 


G2 

G3 

Gi 


-Ih 

Ih 


3 m„ 


+ ..""“"“h' (“Ix' + 12 ™^) 


niq a'lyax AiriqmQajy 


ax 


mQUxx' V 2 


^AA' 


+ Stti^ -\- 


= -Ih 


ax 


+ 
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\/3 V aiv 
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2 
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FA 
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ax 
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_3 
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2 mq Qa\ymqmq 
^ Sm^r 2mcr 5q;| 


(lla^A - 6a"v) 


2mlal 


= -Ih 


'A 

mo 


2 m 


q mq I2mqmq 3 a\ymqmq\ 


3al 


m" 


= -Ih 


2mg 2^mqmQ AmqfriQa^y 

5 «a 
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2. Sturmian Basis 

a. 1/2+ 
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/• 5/2+ 
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APPENDIX D: HADRONIC TENSOR 


The hadronic tensor for these semileptonic decays takes the form 

Hf,^ = -aG^^+P++{p + p')^{p + p')^+P+-{p + p')^{p-p')i, 

+ P-+{P - p')ii(p + P')^ + P — {P- p')ii{p- P')^ 

+ iie-p.^pc,{p + pY{p-pY ■ 

The forms of the terms a, P±± and 7 for the different hnal states we consider are given in the subsections below. 


1 . 1 / 2 + 


a(l/2+) =2{[(toaq - toaJ^ - + [(toaq + toaJ^ - g^]G?} , (Dl) 

1=3 7=3 

/3++(1/2+)= (D2) 

i=l7=l 
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with 

All = A'li = 2 , 

^22 = 2 H'm-AQ + - q'^], 

^'^Aq 

^33 = +”^AJ^ -g^], 

Ag 

Ai2 = (toaq+toa), 

itt-Aq 

A23 = -[(waq+toaJ^-9^], 

■mAQlTlA, 

2 

-431 = -(waq+toaJ, 

mA, 

A 22 = :7^[(’7^AQ-toaJ^-9^], 

^^Aq 

^ -waJ^-9^], 

Ag 

^12 = -^(’T^Aq -WaJ, 
rriAQ 

A 23 = --[(waq-toaJ^-9^], 

rriAQmA, 

2 

7(1/2+) =4FiGi. (D3) 

2. 1/2“ 




45 
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— {mA^+niAj, 
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7 ( 1 / 2 -) =4^^iGi. 


3. 3/2" 


*=4j=4 ^ 

a(3/2-)= ^ 

i=i,j=i 

where Y' = and the non-vanishing coefficients are 

Bii = X[{mAQ+mA^f - q\ 

^44 = 4mi_jTO^J(mAQ+ toaJ^ - g^], 

^14 = = WAQmA,[mi^ - 2 (to^_^ -h g^)TO^Q -h {m\ - q^)% 

B'li = X[{mAQ - ruA^f - q% 

B 44 = 4mi^TO^J(TOAQ -TOAj^-g^], 

*=4j=4 

/3++(3/2-)= ^ _(7l,,F*^/-+4l'4G.G,), 

*=ij=i 

where Y = 12m\^m\^, X = -I- and the Aij are 

^11 = = AXm\^m\^, 

A 22 = Xm^J(mAQ+ toaJ^ - g^], 

A33 = Xml^KmAQ + mA,)"^ - q'^], 

A 44 = 4TOl^miJ(TOAQ + mA,f - g^], 

Ai 2 = AXmAQm\^{mAQ +mAj, 

A 23 = 2XmAQmA,[(mAQ + mA,)^ - g^], 

A 31 = AXm\^mA,{mAQ+mAj, 

Ai4 = -8mi^TOAj(mA<3 +2TOAjg^ + (TOAQ -TOAj(mAQ +toaJ^], 

^24 = 4TOAQmij(TOAQ + mA, f - q^][ml^ - m\^ - q\ 


(D 6 ) 


(D7) 


(D 8 ) 
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^34 = 4migmA,[(TOAQ + mA,)^ - - ml^ - q^], 

A '22 = -^rnlJimAQ - mA,f - q^], 

A '33 = Xml^[(mAQ - niA.f - q^], 

A '44 = 4mi^mlJ(mAQ - mA,f - q^], 

A [2 = 4XmAQmi^(TOA, - toaq), 

A23 = 2XmAQmA^[{mAQ - ttiaJ^ - q"^], 

A'si = 4:Xml^mA,{mA,-mAQ), 

A'44 = 8m\^ml^[{mAQ + 2mAJq'^ - (waq + toaJ(toaq - toaJ^], 

^24 = 4TO^^miJ(mAQ - mA, f - q^][ml^ - m\^ - q\ 

^34 = 4m\^mA,[(mAQ - mA,f - q%m\^ - m\^ - q% 

2 

7(3/2") = —^- —{[{mAQ-mA^f - q\FiGimAQmA^+FiGx[{mAQ+mA,f - q^]) 

+ F 4 G 4 m\^m\^ + ^/GimAgmA, [(toaq + waJ^ - q^]}- 


4. 3/2+ 

i=4,i=4 

a(3/2+)= ^ 

i=l7=l 

where the non-vanishing coefficients are 

Bii = X[{mAQ - niA^f - q% 

B44 = 4 ml^ml^[{mAQ - mA, f - q'^], 

Bi4 = B[4 = mAQmA,[mi^ - 2 {ml^ + q^)ml^ + {m\ - q^)% 
B'li = X[(toaq+toaJ^- g^], 

B'44 = 4 m\^m\JSjnKQ F ruA^f - q% 

i=4:,j=i . 

/3++(3/2+)= ^ _(Al,,.F,i/,-+7l^G,G,), 

*=i.j=i 

where Y = 12toaqTOa^, X = (toa^ -I- TOa, ~ 9^)^ ~ 4771^^441^^, and 

4 I 11 = A ’44 = 4Xmi^mi^, 

A 22 = Xm^J(TOAQ - toaJ^ - g^], 

A 33 = Xml^KmAQ - mA,)"^ - q'^], 

A44 = ‘imX^m'i^limAQ - mA, f - q'^], 

Ai 2 = 4XmAQmA^(TOA, - toaq), 

A23 = 2XmAQmA^[{mAQ - ttiaJ^ - q"^], 

A 31 = AXml^mA^imA^- ttiaq), 


(D9) 


(DIO) 


(Dll) 
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Ai4 = 8 mi^mlJ(mAQ-2mAjq^ - (niAQ+mAj(mAQ 

A 24 = 4ml^mlJ(mAQ - niA.f - - m\^ - q% 

A 34 = 4ml^mA,[{mAQ - mA,f - q^][m\^ - m\^ - q% 

A'22 = XrniJSjnhQ + mA^f - q% 

^33 = Xrni^[{mAQ + ruA^f - q% 

A 44 = 4ml^mlj(mAQ + rriA.f - q\ 

^12 = 4XmAQm\^{mA,+mAQ), 

^23 = 2XmAQmA^[{mAQ + mA^f - q% 

A'^^ = 4:Xml^mA,{mA,+mAQ), 

A[^ = -8ml^ml^[{mAQ - 2mA,)q‘^ - {mAQ - mA,){mAQ - mA,f], 

A '24 = 4m\^m\^[{mAQ + mA,f - q%m\^ - m\^ - q% 

^34 = 4m\^mA^[{mAQ + waJ^ _ q^][nii^ - m\^ - q% 

7(3/2+) = 2^ 2 {[("^Aq + mA, f - q\FiG 4 mAQmA, + J^iGi[(toaq - toaJ^ - 9^]) 

+ F 4 G 4 m\^m\^ + F4GitoaqTOaJ(toaq - waJ^ - 9 ^]}. 

5. 5/2+ 


■i—4,j—4 

a(5/2+)= ^ ±X{B,,F,F,+B[^G.G,), 

*=ij=i ^ 

where the non-vanishing coefficients are 

Bii = X[{mAQ - mA^^ - q\ 

B 44 = 8 m\^m\^[{mAQ + ruA^f - q\ 

Bi4 = B[^ = 2 XmAQmA^, 

B'li = X[{mAQ+mA,Y - q\ 

B 44 = 3mi^TOAj(TOAQ- toaJ^-9^], 

i=4j=4 ^ 

/3++(5/2+)= _(A,,F.F,-+2l^G,G4, 

*=ij=i ^ 

where Yi = 80m® and 

All = A/i =4X'^ml^ml^, 

A 22 = X^mlJ(mAQ+mA,f - q^], 

A 33 = X^ml^KmAQ+mA,y - q^], 

A 44 = 4mAQmX^l(mAQ + iriA^)^ — q^jlq"^ — (2mAQ + mX^)q^ + (iriAQ — m-A^y], 

Ai 2 = 4X^mAQmljmAQ+mA,), 


(D12) 


(D13) 


(D14) 
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A23 = 2 X^mAQmAj(mAQ +mAgf - q^], 

^31 = 4X^mAQmAjmAQ+mA,), 

Ai 4 = -SXmA^mljimAQ+2mAjq^ + (mAQ - ruAjimAQ+mA,f], 
A24 = 4 Xml^mlJ(mAQ + ruA^f - q^][rnA^ - ml^ - q% 

A34 = 4 Xml^mA,[{mAQ + toaJ^ - q‘^][ml^ - m\^ - g^], 

A '22 = X^ml^[{mAQ - mA,f - q% 

A '33 = X^ml^[{mAQ - rriAj^ - q% 

A '44 = 4mi^mlJ(mAQ-mA,f - q^]lq^ - ( 2 ml^+mljq^ + (ml^ 
A'i 2 = 4X^mAQmlJmA, - mAc), 

A 23 = 2 X^mAQmAj(mAQ - niAgf - q^], 

A'31 = 4 X^ 771 ^^ TOA^mA, - mA^), 

A\^ = SXm\^m\J,{mAQ + 2 .mA^)q^ - (toaq +toaJ(toaq - waJ^], 
A 24 = 4Xml^mlj(mAQ - rriA.f - q^][ml^ - m\^ - q% 

As 4 = 4Xml^mAj(mAQ - waJ^ - q^][ml^ - m\^ - q% 


7 ( 5 / 2 +) = 


Aq - 2 wAq(^a, + g^) + {rnl^ + q^f 


^F4G4rnA^m\^ + F4GimAQmA^ (toaq - waJ^ - q^ 


+ [(™Aq + "^aJ^ - q^] {FiG4mAQmA, + FiGi (mAg - toaJ^ - q^ ) | • 


APPENDIX E: CONSTRUCTING HIGHER SPIN REPRESENTATIONS 


It is necessary to construct explicit representations for the spin-3/2 and spin- 5/2 baryons that we treat. In 
the case of the former, the vector-spinor field u‘^{p\ s') must satisfy 


P'aU'^ip'■.s') = 0, 7aW“(p',s') = 0, = 777 ( 3 / 2 )(p', s') • 


(El) 


A suitable representation can be constructed by using the usual Dirac spin-1/2 spinors, together with the 
‘polarization’ vectors e^(p',S 2 ). These vectors satisfy 


pA^iP. Sz) = 0, £*(p', Sz)e+(p , s(,) = -Ss 


(E2) 


Our representation of the spin-3/2 Rarita-Schwinger vector-spinor u^{p', M) is given by the Clebsch-Gordan 


sum 


Ufj,{p',M) = ifj.{p', m)u{p', M — 777)(3/2M|lr77,1/2, M — m). 


(E3) 


This satisfies all of the conditions required. 

A representation of the spin-5/2 spinor u°‘^{p',s) can be constructed in a similar way, but there are two 
additional constraints that must be satisfied. The first is that the spinor must be symmetric in its Lorentz 
indices, and the second is that it must be traceless when the two indices are contracted, i.e. 

<(p',s) = 0. (E4) 

Such a representation can be built in one of two ways. We can use the previously constructed spin-3/2 spinor, 
and the vector e, to write 


. M) E/ib/ 'm)uizip', M - m){5/2M\lm,3/2,M - m). (E5) 


m 
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Alternatively, we can first construct a spin-2 tensor from two of the e vectors as 

,M) = ^ €f_i{p',m)e^{p', M - m){2M\lm, 1,M - m) (E6) 

m 

The symmetry properties of the Clebsch-Gordan coefficients guarantee that this tensor is symmetric in its 
indices. The spin-5/2 spinor is then 

Ufi„{p', M) = ^ A^^{p',m)u{p',M - m){5/2M\2m, 1/2, M - m). (E7) 

m 

These two representations are equivalent, but the manifest symmetry of the second representation allows us to 
see the symmetry in u^i, in an obvious way. 

The conditions 


s') = p'^u°^^(p', s') = 0 (E8) 

are clearly satisfied, since each vector e satisfies p' ■ e=0 (and the second equality also follows from the symmetry 
in the indices). It is easy to check that the auxiliary conditions 

7aM“^(p', s') = 7/3 u“^(p', s') = 0 (E9) 

are satisfied, as are 

s') = m^(5/2)M“^(p', s'). (ElO) 

The traceless condition 

ffa/3M“^(p',s') = 0, (Ell) 

is less obvious, but follows from the tracelessness of Af^^. This, in turn, follows from the symmetry properties 
of the Clebsch-Gordan sum in A^^^, and the properties of the e vectors. 
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